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^ . ABSTRACT 

' We construct aligned and unaligned stationary perturbation configurations in a com- 

QQ I posite system of stellar and coplanarly magnetized gaseous singular isothermal discs 

■ (SIDs) coupled by gravity. This study extends recent analyses on (magnetized) SIDs 
' of Shu et al. , Lou and Lou & Shen. By this model, we intend to provide a conceptual 

framework to gain insights for multi-wavelength large-scale structural observations of 

■ disc galaxies. Both SIDs are approximated to be razor-thin and are in a self-consistent 

axisymmetric background equilibrium with power-law surface mass densities and flat 
rotation curves. The gaseous SID is embedded with a coplanar azimuthal magnetic 

i' field Bg{r) of a radial scaling r~^/^ that is not force- free. In comparison with SID prob- 

lems studied earlier, there exist three possible classes of stationary solutions allowed 
by more dynamic freedoms. To identify physical solutions, we explore parameter space 
involving three dimensionless parameters: ratio A of Alfven speed to sound speed in 
^ , the magnetized gaseous SID, ratio /? for the square of the stellar velocity dispersion 

to the gas sound speed and ratio i5 of the surface mass densities of the two SIDs. 
For both aligned and unaligned spiral cases with azimuthal periodicities |m| > 2, 
one of the three solution branches is always physical, while the other two branches 
^ • might become invalid when /? exceeds certain critical values. For the onset criteria 

I from an axisymmetric equilibrium to aligned secular bar-like instabilities, the corre- 

sponding T/|W — A4\ ratio, which varies with A, P and S, may be considerably lower 
than the oft-quoted value of T/|W| ~ 0.14, where T is the total kinetic energy, W 
is the total gravitational potential energy and A4 is the total magnetic energy. For 
unaligned spiral cases, we examine marginal instabilities for axisymmetric (|m| = 0) 
and non-axisymmetric (|to| > 0) disturbances. The resulting marginal stability curves 
differ from the previous ones. The case of a composite partial MSID system is also 
investigated to include the gravitational effect of an axisymmetric dark matter halo 
on the SID equilibrium. We further examine the phase relationship among the mass 
densities of the two SIDs and azimuthal magnetic field perturbation. Our exact global 
perturbation solutions and critical points are valuable for testing numerical magne- 
tohydrodynamic codes. For galactic applications, our model analysis contains more 
realistic elements and offer useful insights for structures and dynamics of disc galaxies 
consisting of stars and magnetized gas. 
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1 INTRODUCTION 

In galactic contexts, we venture to formulate a theoretical 
magnetohydrodynamic (MHD) disc problem to explore pos- 



sible large-scale structures and dynamics of stationary MHD 
density waves in a composite system of stellar and magne- 
tized interstellar medium (ISM) gas discs. The two grav- 
itationally coupled discs are treated as 'fluid' and 'mag- 
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netofluid' respectively and are both expediently approxi- 
mated as razor-thin singular isothermal discs (SIDs) with 
the gaseous SID being embedded with a coplanar azimuthal 
magnetic field. For the gravitational eflFect of a massive ax- 
isymmetric dark matter halo, we prescribe a background 
composite system of two coupled partial SIDs (Syer & 
Tremaine 1996; Shu et al. 2000; Lou 2002; Lou & Shen 2003; 
Shen & Lou 2003). In our model analysis, we construct sta- 
tionary aligned and unaligned logarithmic spiral MHD per- 
turbation configurations in a composite system of two SIDs 
with flat rotation curves, and attempt to relate various mor- 
phologies of disc galaxies, including barred and lopsided, 
barred and normal spiral structures. For possible observa- 
tional diagnostics, we derive phase relationships among per- 
turbation patterns of the stellar surface mass density, the 
gas surface mass density and the azimuthal magnetic field. 

This introduction serves two purposes. The first one is 
to provide a general background information relevant to the 
problem at hand and the second one is to give reasons of 
pursuing this MHD disc problem. 

In a pioneering study of a composite system of stellar 
and gas discs coupled by gravity, Lin & Shu (1966, 1968) 
used a stellar distribution function and a gas fluid disc de- 
scription to derive and analyze the local dispersion relation 
of galactic spiral density waves. Since then, there have been 
extensive theoretical studies on perturbation configurations 
and stability properties of a composite disc system, mainly 
in galactic contexts. Kato (1972) investigated oscillations 
and overstabilities of density waves using a formalism sim- 
ilar to that of Lin & Shu (1966, 1968). In a two-fluid for- 
malism. Jog & Solomon (1984a, b) examined the growth of 
local axisymmetric perturbations in a composite disc sys- 
tem. Bertin & Romeo (1988) studied the influence of a gas 
disc on spiral modes in a two-fiuid model framework. Van- 
dervoort (1991a, b) studied the influence of interstellar gas 
on oscillations and stabilities of spheroidal galaxies. The 
two-fluid approach was also adopted in a stability study of 
a two-component disc system with flnite disc thickness by 
Romeo (1992). The analysis for morphologies of disc galax- 
ies was performed by Lowe et al. (1994). For the stability of 
a composite disc system, different effective Qef / parameters 
(Safronov 1960; Toomre 1964) have been suggested using a 
two-fluid formalism by Elmegreen (1995) and Jog (1996). 
Recently, Lou & Fan (1998b) used the two-fluid formalism 
to study properties of open and tight-winding spiral density- 
wave modes in a composite disc system. Lou & Shen (2003) 
studied stationary global perturbation structures in a two- 
fluid system of SIDs and, instead of a redeflnition of a dif- 
ferent Qeff parameter, Shen & Lou (2003) offered a more 
practical Z)— criterion for the axisymmetric instability in a 
composite SID system. 

A rich class of disc problems involves stability properties 
of SIDs. There have been numerous studies on this subject 
since the pioneering work of Mestol (1963) (e.g. Zang 1976; 
Toomre 1977; Lemos, Kalnajs & Lynden-Bell 1991; Lynden- 
Bell & Lemos 1999; Goodman & Evans 1999; Charkrabarti, 
Laughlin & Shu 2003). Specifically, Syer & Tremaine (1996) 



made an important breakthrough to derive semi-analytic so- 
lutions for stationary perturbation configurations in a class 
of SIDs. Shu et al. (2000) obtained stationary solutions for 
perturbation configurations in an isopedically magnetized 
SID with a flat rotation curve. Through numerical explo- 
rations, they interpreted these stationary aligned and un- 
aligned logarithmic spiral configurations as onsets of bar- 
type and barred-spiral instabilities (see also Galli et al. 
2001). Different from yet complementary to the analysis of 
Shu et al. (2000), Lou (2002) performed a coplanar MHD 
perturbation analysis in a single background SID embedded 
with an azimuthal magnetic field, from the perspective of 
stationary fast and slow MHD density waves (FMDWs and 
SMDWs; Fan & Lou 1996; Lou & Fan 1998a). Lou (2002) 
also derived a form of magnetic virial theorem for an MSID 
and suggested the ratio of rotation energy to the sum of grav- 
itational and magnetic energies to be crucial for the onset of 
bar-like instability in an MSID system. In galactic contexts, 
it would be more realistic to consider large-scale structures 
and dynamics in a composite system of stellar and magne- 
tized ISM discs. As a first step, Lou & Shen (2003) made a 
foray on this model problem, constructed stationary aligned 
and unaligned logarithmic spiral configurations in such a 
composite SID system and further examined axisymmetric 
instability properties (Shen & Lou 2003). 

In disc galaxies, the ISM disc is magnetized with the 
magnetic energy density being comparable to the energy 
densities of thermal gas and of relativistic cosmic-ray gas 
(e.g. Lou & Fan 2003). Information of galactic magnetic 
fields can be estimated by synchrotron radio emissions from 
spiral galaxies. For such a magnetized composite system, 
MHD will play an indispensable role and reveal more real- 
istic aspects of dynamic and diagnostic information. These 
important problems (Shu et al. 2000; Lou 2002; Lou & Shen 
2003) are not only interesting by themselves, but also serve 
as necessary steps for establishing an even more realistic 
model. Motivated by this prospect (Lou & Fan 1998b; Lou 
2002; Lou & Shen 2003), we construct here stationary per- 
turbation configurations for aligned and unaligned logarith- 
mic spiral cases in a composite system of a stellar SID and 
a coplanar ly magnetized gaseous MSID, and discuss their 
stability properties. 

We adopt a relatively simple formalism for a compos- 
ite system involving fiuid and magnetofluid discs coupled by 
gravity. We provide in Section 2 an MHD description for the 
coplanarly magnetized gaseous MSID, obtain conditions for 
background axisymmetric equilibrium state for both stellar 
SID and gaseous MSID, and derive linearized equations for 
coplanar perturbations. There exist aligned and unaligned 
classes of global MHD perturbation solutions; they are ana- 
lyzed in Section 3 and Section 4, respectively. The exact so- 
lutions of stationary perturbations, their stability properties 
and their corresponding phase relationships among pertur- 
bation variables are examined and summarized in Section 5. 
Details are included in Appendices A— E. 



Configurations of composite MSID 3 



2 FLUID-MAGNETOFLUID FORMALISM 

It would be physically more precise to adopt a distribution 
function formalism in dealing with a stellar disc especially in 
terms of singularities and resonances (e.g. Lin & Shu 1966, 
1968; Binney & Tremaine 1987). For the present purpose of 
modelling large-scale stationary perturbation structures and 
for mathematical simplicity (Lou & Shen 2003), it suffices 
to start with the fluid-magnetofluid formalism, including an 
MHD treatment for the gaseous MSID (Lou 2002; Lou & 
Fan 2003). In this section, we present the basic equations 
for the fluid-magnetofiuid system consisting of a stellar SID 
and a gaseous MSID. For flat rotation curves, conditions 
on the background rotational equilibrium with axisymme- 
try can be derived. We then obtain the linearized equations 
for coplanar MHD perturbations in the composite MSID 
system. 

2.1 Basic Nonlinear MHD Equations 

The two SIDs, located at z — 0, are both approximated 
as infinitesimally thin. In our fluid-magnetofluid treatment, 
the two SIDs are coupled through mutual gravitational in- 
teraction. For large-scale stationary perturbations, diffusive 
processes such as viscosity, ambipolar diffusion and ther- 
mal diffusion etc. are ignored. For physical variables under 
consideration, we shall use superscript or subscript s to in- 
dicate an association with the stellar SID and superscript 
or subscript g to indicate an association with the gaseous 
MSID. In cylindrical coordinates (r, 6, z), the basic fluid- 
magnetofluid equations for a composite MSID system can 
be readily written out. 

In the fluid approximation for a stellar SID, the mass 
conservation, the radial component of the momentum equa- 
tion and the azimuthal component of the momentum equa- 
tion are given below in order, namely 

as- 1 a(rs-»-) 1 9(s-j-) _ 

dt ^ r dr 86 ' ^ ' 
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where u" is the radial component of the stellar bulk velocity, 
= rv" is the stellar specific angular momentum along 
the z direction; and v" is the azimuthal component of the 
stellar bulk velocity; </> is the total gravitational potential, 11° 
is the vertically integrated (effective) pressure (sometimes 
referred to as the two-dimensional pressure), and is the 
vertically integrated stellar mass density (i.e. stellar surface 
mass density). 

In the magnetofiuid approximation for the gaseous 
MSID, the mass conservation, the radial component of the 
momentum equation and the azimuthal component of the 
momentum equation are given below in order, namely 
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where is the radial component of the gas bulk velocity, 
f = rv^ is the gas specific angular momentum along the 
z direction, is the azimuthal component of the gas bulk 
velocity, 11^ is the vertically integrated gas pressure (some- 
times referred to as the two-dimensional gas pressure), E^ 
is the vertically integrated gas mass density (i.e. gas surface 
mass density) and Br and Bg are the radial and azimuthal 
components of magnetic field B. The last two terms on the 
right-hand sides of equations ^ and Q are the radial and 
azimuthal components of the Lorentz force due to the copla- 
nar magnetic field. The coupling of the two sets of fluid and 
magnetofiuid equations (1) — (3) and (4) — (6) is effected 
by the total gravitational potential (f> through the Poisson 
integral, namely 

-G(Ss + E^)CrfC 



F(f,{r,d,t) = <p dill 



(7) 



/ [C +r'^ - 2(r cos(7/) - 6)] 
where F is a constant ratio with < F < 1; F(f> is the grav- 
itational potential from the stellar SID and gaseous MSID 
together, and the fraction {1 — F)(j) is attributed to an ax- 
isymmetric dark matter halo that is unresponsive to copla- 
nar MHD perturbations in the composite MSID system (e.g. 
Shu et al. 2000). 

The divergence-free condition for the coplanar magnetic 
field B = {Br,Be, 0) is 
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Or 
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and the radial and azimuthal components of magnetic in- 
duction equation are 
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Equations (1) — (10) form the basis of our theoretical anal 
ysis. 



2.2 Rotational MSID Equilibrium 

For a fiat rotation curve in a stellar SID, we write the cor- 
responding angular rotation rate in the form of 



Q.s(r) = as 



(11) 



where Os, mimicking an 'effective isothermal sound speed', 
represents the velocity dispersion of the stellar SID and Da 
is a dimensionless parameter for stellar SID rotation. We 
invoke the expedient polytropic approximation 

Ho = (12) 
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to relate the two-dimensional pressure and the surface mass 
density. The epicyclic frequency Ka of the stellar disk is de- 
fined by 



, — » . (13) 

r dr 

In parallel, we have the angular rotation rate of the gaseous 
MSID as 

n,{r)=a,^, (14) 

where ag is the isothermal sound speed of the gaseous MSID 
and Dg is a dimensionless parameter for MSID rotation. 
The polytropic relation and the definition of the epicyclic 
frequency Kg for the gaseous MSID are simply 

2^9 



Kg — 2D,g 



(15) 



(16) 



To avoid the magnetic field winding dilemma in a rotat- 
ing disc (e.g. Lou & Fan 1998a), the background coplanar 
magnetic field, which is not force-free, is taken to be purely 
azimuthal about the symmetry axis; 



Bg{r) = 



(17) 



where is a constant (Lou 2002) proportional to the encir- 
cled magnetic flux within r, and 



Br = B^=0 



(18) 



In galactic model applications, one needs to invoke a cen- 
tral bulge or other processes to avoid the divergence of 
Be as r 0. From the radial momemtrum equations 
lj2J and in a rotational equilibrium and the fact that 
Fd(j)/dr = 27rG(Eg + Eg) by Poisson integral Q, one de- 
rives the following expressions for the background surface 
mass densities, namely 

,ai{l + Dl) 1 



2TvGr (1 -I- S) 
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(20) 
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where 5 = Eq/Eq is the surface mass density ratio of the 
two coupled background SIDs and Ca is the Alfven wave 
speed in the MSID defined by 



(21) 



(22) 



dzBe/(47rEo) . 

From equations 1191 and I20I I. it then follows that 

al{l + Dl)=al{l + Dl)-C\/2 . 

Physically, condition 1221 results from the basic fact that 
the same total gravitational force d(j)/dr, including the con- 
tribution from the dark matter halo, acts on both the stellar 
SID and magnetized gaseous SID, and is very useful in our 
analysis below. It should be noted that the rotation rates 
of the two SIDs are different in general (Lou & Shen 2003; 
Shen & Lou 2003) . In dimensionless form, condition 1221 can 
be written in the form of either 



^ n 



/3(1 + D') - 1 + A72 
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D! = ^{l + Dl-Xy2)-1 
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where parameter P = a^/al stands for the square of the 
ratio of the stellar velocity dispersion to the sound speed of 
the MSID, and parameter A'^ = C\/a?g stands for the square 
of the ratio of the Alfven speed to the sound speed in the 
MSID. In disc galaxies, the stellar velocity dispersion is 
usually higher than the sound speed ag, we naturally focus 
on the case of /? > 1 (e.g. Jog & Solomon 1984a, b; Bertin & 
Romeo 1988; Jog 1996; Elmegreen 1995; Lou & Fan 1998b; 
Lou & Shen 2003; Shen & Lou 2003). In this /? > 1 regime, 
it follows from condition 1241 that 

1 + dI < 1 + dI~ \^/2 < 1 + dI , (25) 

implying that for 13 >1, 

d! < Dl . (26) 

Inequality 12611 is very important to identify physically valid 
mathematical solutions of for stationary MHD perturba- 
tions. For specified parameters /3 and A'^, and are re- 
lated to each other linearly by condition 12311 . We emphasize 
that mathematical solutions of Dg and become unphysi- 
cal for either < or Dj. < or both. The key here is that 
by inequality 1261 , we only need to consider Di > because 
Dg must also be positive. In our analysis, we mainly use 
equation 1231 to derive a cubic algebraic equation in terms 
of Da and examine solution properties. 



2.3 Perturbations in a Composite MSID System 

For small coplanar MHD perturbations in a composite MSID 
system, basic nonlinear equations ITl — 1101 can be linearized 
in a straightforward manner, namely 
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for coplanar hydrodynamic perturbations in a stellar SID, 
and 
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for coplanar MHD perturbations in a gaseous MSID, 



F(t>i 



G <p dip 
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for the linearized Poisson integral, and 
d{rhr) dbg 



11/2 
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1 d 



(33) 



(34) 



(35) 



(36) 



for the linearized divergence-free condition and the lin- 
earized magnetic induction equation. 

We do not consider vertical variations along z direc- 
tion across the composite MSID system. With a harmonic 
exp{iujt — imO) dependence for all perturbation variables, 
we introduce complex radial variations ^jLs{t), iJ.g{r), Us{r), 
Ugir), J,(r-), Jg{r), V(r), R{r) and Z{r) for Ef, ul, 
'^it ill ill 01 1 and be, respectively. Thus, hydrodynamic 
equations 1271 — 1291 can be reduced to the form of 

i(Lu - mQs)f^s + --^{rE'oUs) - '-^Js = , (37) 
r or ' 



3 ALIGNED MHD CONFIGURATIONS 

Coplanar perturbations in a composite MSID system can be 
classified as 'aligned' and 'unaligned' solutions (e.g. Kalnajs 
1973; Shu et al. 2000; Lou 2002; Lou & Shen 2003). For 
aligned configurations, all streamlines and magnetic field 
lines are aligned in a composite MSID system. For unaligned 
spiral configurations, neighbouring streamlines shift relative 
to each other in a systematic manner (Kalnajs 1973); the 
same physical scenario holds true for neighbouring magnetic 
field lines (Lou & Fan 1998a). 

In this section, we obtain the stationary dispersion re- 
lation for aligned coplanar MHD perturbations (both full 
and partial MSIDs) by perturbation equations (I33II and 
H37|l — H45|l in the preceding section. The solution behaviours 
and the corresponding phase relationships are analyzed, 
mainly in the context of a full SID system (i.e. F = 1). At 
the end of this section, we derive the MHD virial theorem 
for a composite MSID system and suggest the onset crite- 
rion for secular bar-like instabilities in a composite MSID 
system (Ostriker & Peebles 1973; Binney & Tremaine 1987; 
Shu et al. 2000; Lou 2002). 

3.1 Dispersion Relation for Aligned Perturbations 
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for the stellar SID, where <I>s = a^^s/Eo + V. Similarly, 
MHD equations CT-lCT can be cast into the form of 
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for coplanar MHD perturbations in the gaseous MSID. By 
setting angular frequency ci; = in coplanar MHD perturba- 
tion equations l|33^ and l|37|l — H45|l . we can construct global 
stationary MHD perturbation configurations in a compos- 
ite system of MSIDs without invoking the WKBJ or tight- 
winding approximation and analyze their properties. 



To construct stationary perturbation configurations of 
MSID that are aligned, we set a; = in equations 13711 — 14511 . 
Let us first set d / dt = G ox uo = in equations 1401 — 1451 1 
for coplanar MHD perturbations in the MSID to obtain 
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For aligned perturbations, we take the following potential- 
density pair (Shu et al. 2000; Lou 2002; Lou & Shen 2003) 



such that 
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For a constant iUg, as will be shown presently, combinations 
of equations 11461 and 14811 and expressions 15111 — 1531 1 give 



Jo 



Y9 

^0 



which gives the stationary dispersion relation for aligned 
coplanar perturbations in a stellar SID alone, that is, a single 
SID without magnetic field. Substituting expressions IIH 
and Ijl9|l of Q,s and Eq (with F = 1) into equation H68|l . we 
obtain 



C7i/(2n,r)-n,r 



(56) 



As $g is constant by equation I54L it is clear that iUg is an- 
other constant. Consequently, equations 15511 and 1471 give 
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Substitutions of expressions 1531 , 15511 and 1561 into equation 
1571 give a relation between yLg and namely 
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where two coefficients K. and A are defined by 
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respectively. In parallel, we set lj = in equations l|37^ — Ij39|l 
for coplanar perturbations in the stellar disc and use rela- 
tions 15111 — 1531 1 to obtain* another relation between fis and 
Hg, namely 
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Combining equations 16111 and 1581 , we derive the stationary 
dispersion relation 
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(62) 

for aligned coplanar MHD perturbations in a composite 
MSID system. In the absence of the mutual gravitational 
coupling between the MSID and the stellar SID, represented 
by the term on the right-hand side of equation 1621 , the left- 
hand side of equation l|62|l would give rise to two separate 
dispersion relations, one for the stellar SID and one for the 
gaseous MSID. The first one would be 



(m^ - 2)n'ir^ + 2nG\m\T.or - 



2 2 n 
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* In view of the exchange symmetry between the stellar and gas 
SIDs, relation ED 

can also be obtained by simply setting Ca = 
in equation I58i and switching subscripts g and s. 
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which is simply equation (26) of Shu et al. (2000). 

The second factor in the curly braces on the left-hand 
side of equation 16211 is 
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(65) 

which is the stationary dispersion relation for aligned copla- 
nar MHD perturbations in the gaseous MSID alone, that is, 
a single MSID with a coplanar magnetic field. Equation l|65^ 
is simply equation (3.2.6) of Lou (2002). 

By the above results, it is clear that equation 16211 rep- 
resents the dispersion relation for stationary MHD density 
waves in a composite MSID system. As both SID and MSID 
rotate, the stationarity of MHD density wave patterns in an 
inertial frame of reference imposes conditions on dimension- 
less rotation parameter (or equivalently, -Dg). For the 
following analysis, we substitute expressions lllH . 114L 1191 1 
and I2OI 1 into equation 1621 to yield another form of station- 
ary dispersion relation in a composite full or partial MSID 
system, namely 

1 1 



■m?Dl 



{rn -2)DI + F\m\ 



- A2(m2 

l + Dl _ 
1+5 



1/2) 



4 r-,4 

m Vg 



2Dl + 



A-' + l-F- 



l + Dl- \'/2 5 



+ 



2 \ 2 
xm A 



l+S 
AV2 5 



2\' 



\m\ l + <5 



X raDi 



A^ 

4 



AV2)5 



(1 + 5)2 



(66) 

A substitution of expression ^ for Dl into ^ gives a 
cubic algebraic equation in terms of D\. Or equivalently, a 
substitution of expression for D1 into would yield 
a cubic algebraic equation in terms of Dg. As noted earlier, 
we focus on the cubic equation of D\ to identify physical 
solutions of D1 > 0, because < Dg as a result of > 
in typical disc galaxies. 

For a later examination of spatial phase relationship 
between azimuthal magnetic field perturbation be and the 
surface mass density perturbation Ef of the gaseous MSID, 
we combine equations l|49^ and H50^ to obtain 
^ _ iUgBs 



2mQ,gr 



(67) 



relating the ^-component of the magnetic field perturbation 
Z and the radial gas flow speed perturbation iUg. Using 
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equations 1571 and 1671 1. one can eliminate iUg to obtain 

^_ Be^,, {2Q.y-C\/2) 

2Eg [m2n2r2 -C2(m2 - 1/2)] ' ^ ' 

that can be further reduced to 

_ 2Sg [m^Dl - X^m^ - 1/2)] 
^ Be 2D2_A2/2 ■ ^ 

For real fJ^g/Z, the sign of the right-hand side of the equation 
16911 will determine the phase relationship between the az- 
imuthal magnetic field perturbation bg and the gas surface 
mass density perturbation in the gaseous MSID. That is, 
bg and Ej are in and out of phase for a positive and negative 
right-hand side of the equation 1691 , respectively. 

3.2 Axisymmetric Disturbances with \m\ = 

For aligned axisymmetric disturbances with \m\ = 0, it 
would be inappropriate to directly use relation 16211 or I66II . 
One should carefully examine equations H37|l — l|45|l with 
Lu = m ^ 0. With Ug ^ Us ^ R ^ 0, equations I^ZJ, 
13911 . 14UII . 14211 . 1431 1 and 14511 can be satisfied, and equation 
14411 is identically zero. By choosing Z oc r^^^^, jig oc r^^, 
^3 oc r~^, Jg <yi r, Js r, V <yi Inr, ^g oc Inr-f constant 
and oc Inr-f- constant, the two remaining equations 13811 
and 1411 are consistent with a rescaling of the axisymmetric 
background. In the present context, this rescaling is some- 
what trivial. We turn to cases of |m| > 1 below. 

3.3 Nonaxisymmetric Disturbances with |m| > 1 

3. 3. 1 Behaviours of solutions 

We first come to the aligned case of jmj = 1. For a full 
composite MSID system with _F = 1, it is easy to verify that 
equation 16611 can be satisfied for arbitrary Dl. This is quite 
similar to cases of a single SID studied by Shu et al. (2000) , 
of a single MSID studied by Lou (2002) and a composite SID 
system studied by Lou & Shen (2003). However, we note that 
for a single partial MSID (Lou 2002) and for a composite 
system of two coupled partial SIDs (Lou & Shen 2003), such 
stationary aligned eccentric |mj = 1 perturbations are not 
allowed for arbitrary . Likewise, in our case of a composite 
partial MSID system with F < 1, it is easy to see that 
equation 16611 can no longer be satisfied for arbitrary by 
simply setting \m\ — 1. 

In the following, we analyze cases of \m\ > 2 and focus 
on the special case of _F = 1 for a full composite MSID sys- 
tem. As mentioned earlier, we work in the parameter regime 
of /3 > 1 with a typical disc galaxy in mind. By an exten- 
sive numerical exploration of aligned cases from \m\ = 2 to 
\m\ — 5, we note empirically that cases of jml > 2 are fairly 
similar to the case of |m| = 2. For this reason, we shall 
mainly consider the case of \m\ = 2. As noted earlier, we 
can substitute expression 12311 for Dg into equation 16611 to 
derive a cubic algebraic equation of y = as 

A + By + Cy^ +Vy^ ^0 , (70) 



m|=2, 5=0.2, 1^=0.09 



1.5- 




1 1.5 2 2.5 o 3 3.5 4 4.5 5 



Figure 1. Three solution curves of versus /3 for the aligned 
case with |m| = 2, 5 = 0.2 and A2 = 0.09. 

where jmj = 2 and F = 1, and the four coefficients A, B, C 
and D are explicitly defined by 

^ = 24 - 32 /? ~ 3 a" + 12 -I- 48 (5 - 80 /3 (5 
-6 P X^S + 8 0^ + 6 + 32 - 6 X'^S , 

B = -12 - 3 /? - 48 + 32 /3 -f 48 
-40 13 5 - 245 + 6 X"^ - 6X^5 + 3X^5 , 

C = 3 P X'^5 + 40 P 5 + 64 f3 - 24 13^ 
and 

V = -16/3^5-16/3^ , 
respectively. 

Formally, three mathematical solutions for y = can 
be written out analytically from cubic equation 1701 or from 
the more general cubic dispersion relation 1661 with jmj > 2 
and F < 1, although the solution expressions are fairly in- 
volved (see Appendix C for details). Practically, we explore 
numerically various parameter regimes for the three DI so- 
lutions. For specified values of parameters 5 and A^, we show 
three mathematical solutions y = Dl versus /? in Figs.0-|3 

For typical parameters, the three solution branches of 
do not intersect with each other. For the convenience of 
discussion, we use j/i , j/2 and j/3 to denote the upper, middle 
and lower solution branches, respectively. Generally speak- 
ing, the upper two solution branches, t/i and 1/2, are quali- 
tatively similar to those in a composite unmagnetized SID 
system (Lou & Shen 2003), that is, j/i branch remains always 
positive, while t/2 decreases monotonically with increasing /3 
and becomes negative when P exceeds a certain critical value 
Pel- As expected for small A'^ (i.e. weak magnetic field; see 
Fig. , j/i and 1/2 respectively obey the limits of those in a 
composite SIDs without magnetic field [see equations (54) 
and (55) of Lou & Shen (2003)], showing that as A^ 0, yi 
and j/2 consistently approach the upper and lower branches 
respectively of Lou & Shen (2003). One novel feature is the 
lowest solution branch 3/3 owing to the presence of magnetic 
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in 

Q 0.5 





|m|=2, 5=0.2, ?i^=1 










Pc2 ^ 


^^"~yr~~- — — - — — ; 



1 1.5 2 2.5 p 3 3.5 4 4.5 5 

Figure 2. Three solution curves of versus /3 for the aligned 
case with \m\ = 2, 5 = 0.2 and = 1. 




Figure 3. Three solution curves of versus /3 for the aligned 
case with \m\ = 2, 5 = 0.2 and = 3.61. 



field. In most cases, j/3 is negative when /3 varies from 1 
to +00, but for some special parameters, j/3 may become 
positive when /3 becomes smaller than a critical value /3c2. 
Analytical expressions of (id and /3c2 for the case of \m\ = 2 
can be derived from a quadratic equation and are given by 
the following pair 
1 



/3ci = 



[40 5 + 16 + 3 



and 

/3c2 



8 + 325 

(64 5^ + 128 5 - 144 X^S'^ + 64 - 192 
+201 A*(5^ + 24 A'' + 144 A*5 - 48 X^f^'^] 



:[405 + 16 + 3A^5 



(71) 



8 + 32(5 
(64 5' 

+201A''5' + 24 A* + 144 A''^ 



128 5 - 144 A'5' + 64 - 192 A'5 



(72) 



respectively, where fid remains always larger than /3c2. We 
note the following. First, there is no essential mathematical 
difficulty of obtaining more general forms of fid and /3c2 for 



arbitrary |m| values with < F < 1. Secondly, by setting 
a' = in expression I7H for fid, we obtain 

1 



1 + 



(73) 



45 + 1^ 

consistent with the case of |mj = 2 in expression (56) for 
fic by Lou & Shen (2003) . Finally, while /3c2 is real in many 
cases [i.e. a positive determinant in both expressions 1711 
and Ij72|l ]. /3c2 may be too small to be discernible in the pa- 
rameter regime of /3 > 1. For some special parameters spec- 
ified, Pc2 becomes noticeable. For example, in Fig. |21 with 
5 — 0.2 and A' — 1, we have /3c2 = 1.1667 > 1 according to 
expression In this case, for /3 smaller than 1.1667 (still 
quite restrictive), there exist three positive solutions of D^. 
For Pc2 < fi < /3ci, the upper two branches j/i and y2 are 
positive, corresponding to two possible stationary perturba- 
tion modes. When /3 exceeds Pd, only yi remains positive, 
corresponding to one possible stationary perturbation mode. 
For the usual case of /3c2 < 1, there are at most two possible 
stationary modes when /3 ^ 1 (see Fig.lSJ. 

3.3.2 Phase relationships among perturbation variables 

We now examine phase relationships among the azimuthal 
magnetic field and the surface mass density perturbations, 
because they may provide clues for magnetized spiral galax- 
ies through optical and synchrotron radio observations (e.g. 
Mathewson et al. 1972; Beck & Hoernes 1996; Fan & Lou 
1996; Lou & Fan 1998a, 2002, 2003; Frick et al. 2000, 2001; 
Lou 2002; Lou et al. 2002). For the phase relationship be- 
tween the two surface mass density perturbations and 
/is, a combination of expressions lllll . (1191 and equation 
1611 gives 

l^, _ 1 [y(m'-2)-m'](l + 5) 

Ms l"i|(j/ + l) 

where y = D^. A substitution of expression 1741 into equa- 
tion 1701 leads to a cubic algebraic equation of /ig//is. We 
show different curves of /-ig/ns versus /3 in Figs. HI and |^ by 
specifying different values of 5 and A^. 

In parallel, we examine the phase relationship between 
the azimuthal magnetic field perturbation be and the sur- 
face mass density perturbation in the gaseous MSID. By 
equation 1691 . we introduce a dimensionless q parameter 



[m^Dl -\'{m' -1/2)] 



2DI 



AV2 



(75) 



whose sign determines the phase relationship between Z and 
fig. Using equation 1751 . we may express Dg as a function 
of q, namely 



A' (g + 2m' - l) 



4g + 2m2 ' ^^^^^ 
A combination of equation 1661 and expression 1241 gives 
a cubic algebraic equation of D^; and a substitution of ex- 
pression 1)76^ into the resulting cubic equation of leads 
to a cubic equation in terms of q (with jmj = 2). We present 
different curves of q versus fi for different values of 5 and A' 

in Figs.inmi 
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Figure 5. Curves of ^lg/^ls versus /3 for the aligned \m\ = 2 case Figure 7. Curves of fig/Z versus /3 for the aligned |m| = 2 case 

with 5 = 0.2 and = 3.61. The uppermost branch related to yj, ^jt^ S = 0.2 and = 3.61. The uppermost branch related to ys 

is only discernible for /3 < 1 and is not shown here. Specifically, appears somewhat fiat at about 44 and is not shown here, 
this is caused by a negative j/3 when /3 > 1 (see Fig.E}. 



From expressions 1741 and 1751 , one can show that both 
larger /ig//is and fig/Z correspond to smaller y. The lowest 
branches of jj,g/fis or i-ig/Z are then related to the uppermost 
j/i; the middle branches of fJ^g/fJ-s or fJ.g/Z are related to the 
middle 1/2; and the uppermost branches of i-ig/i-is or i-ig/Z 
are related to the lowest j/3 as shown in Figs. |1J-(7| 

Several interesting features are noted here for these 
curves. First, when is small (see Fig. 2J, the lowest and 
middle branches of fig/ Us respectively follow those in a com- 
posite SID system without magnetic fields [see equations 
(58) and (59) of Lou & Shen (2003)], showing consistently 
that as 0, the middle and lowest branches of fig/fJ-s 

correspond to the lower and upper branches of Lou & Shen 
(2003), respectively. 

Secondly, fj.g/fis and fJ^g/Z of each solution branch of y 
bear definite signs in most cases. For the case of 5 = 0.2 
and A^ = 0.09 as shown in Figs. |1| and |S| the stationary 
perturbation mode of yi branch has a phase relationship 



between the surface mass density of the two SIDs (i.e. /ig//ia) 
and a phase relationship between the surface mass density 
of the MSID and the azimuthal magnetic field (i.e. Hg/Z) 
being both always out of phase. In most cases, the stationary 
perturbation mode of 7/2 branch has a fig/fis being in phase 
(see Fig.QJ and a, jj,g/Z being out of phase (see FigEJ. The 
stationary perturbation mode of yj, branch has jJ.g/i-ia and 
l^g/Z being both always in phase (see Figs.0]andO. 

Thirdly, for the stationary perturbation mode of 2/2 
branch, both fig/i-ia and fig/Z have a common zero point 
at Pco that increases with increasing A^. For a sufficiently 
large A^ (e.g. A'^ = 3.61), /3cO becomes greater than 1. Thus, 
fj-g/fj-s and ^J.g/Z may become zero or even carry opposite 
signs (see Figs. |K| and |7|l. Our analysis show that the middle 
2/2 branch for fJ-g/ns and fig/Z has a common root at 

(6 + 3(5)A2 + 8<5+ 16 



(77) 



24(1 + ^) 

For a fixed S, /3cO increases with increasing A^. Specifically 
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for 5 = 0.2 and A'^ = 0.09 in expression 1771 . we have 
PcO = 0.6317 smaller than 1. Therefore, the middle 1/2 branch 
remains always positive (see Fig.|lJ. On the other hand, for 
= 3.61 in expression l)77|l . we have /3co = 1.4384 larger 
than 1 (see Fig.l^. 

It is not surprising that the middle y2 branch of fig/^a 
and ^ig/Z has a common root as result of /ij = 0. The special 
Hg — case may be a possible situation in our model. By 
equation 1551 1. it follows that Jg vanishes. As fis and Z do 
not vanish, we see that V, $g and consequently Ug do not 
vanish by expressions 15311 . 1541 and 1561 . 

Finally, our analysis show that the uppermost branch 
of Mff/^s and iJ.g/Z corresponding to the 1/3 branch has a 
root at /3 = 0. It follows that j/3 remains always positive. 



3.4 Secular Barlike Instabilities 

For the onset of aligned secular barlike instabilities in a com- 
posite MSID system, we here derive the MHD virial theorem 
for a composite MSID system (either full or partial) from the 
background radial equilibrium conditions Q and Q 



W = 



r(Eo + Eg) — 27rrdr = R, (87) 



M 



2r 



2nr dr ■ 



2G(i{l + 5) 



R 



(88) 



(78) 



and 



Adding equations l)78|l and Ij79|l . we readily obtain 



(ESn^ + Egn^)r = £-(a?E5+a^Eg) + (ES+Eg)^ + Eg^. 

(80) 

Multiplying equation 1801 by 2-Kr^dr and integrating from 
to a finite radius R, that is allowed to approach infinity 
eventually, we have the following MHD virial theorem. 

2{T + U)+W-M= 2iTR'[alT,l{R) + alT.l{R)] (81) 

within radius R, where 



Based on early numerical simulations (Hohl 1971; Miller et 
al. 1970), it was known that a thin self-gravitating disc in 
rotation may rapidly evolve into bar configurations (e.g. Bin- 
ney & Tremaine 1987). Ostriker & Peebles (1973) suggested 
an approximate criterion T/|W| < 0.14±0.02, necessary but 
not sufficient, for stability against bar- type instabilities, on 
the basis of their A'^— body numerical explorations involving 
300 particles. In the presence of a coplanar azimuthal mag- 
netic field with a radial scaling of oc r~^/^, Lou (2002) pro- 
posed that the ratio T /\W — M\ may play the role of T 
in an unmagnetized SID to determine the onset of instabil- 
ity in a single MSID. By the above analogy and a natural 
extension, we use T /\'W — M\ instead of T/|VV| to examine 
the onset criterion of instability in a full composite MSID 
system in the presence of coplanar non-axisymmetric aligned 
MHD perturbations. Using expressions 1861 . 1871 and 1881 . 
we arrive at 

T _ -2S + 25f3 + 5\^ + {25p + 2l3)Dl 
\'W -M\ ^ 25A2 -1-4/3 + 45/? + (4 /3-f45/3)D2 

2 5A2 + 2/3(1 + 5)(Lif + l) ' ^ ' 

indicating that the value of T/|VV — M\ falls between 
and 1/2 as in the unmagnetized case (Lou & Shen 2003) 
and increases with increasing Dl. Therefore in a composite 
MSID system, the three possible values of correspond to 
three different values of T/|VV — M\ ratio; and larger values 
of correspond to higher T/jW — M\ ratios. 

We here examine a few cases to illustrate the utility 
of our proposed onset criterion for instability. First, we set 



T= r -T.l{rnsf2T^Tdr + -E'o{rQgf2Trrdr (82) i'"! = 2, 5 = 0.2, /3 = 1 and X'' = l. From equation COJ, 
Jo "2 Jo '2 one obtains three D? solutions D? = 1.8246, D? = 0.5513 



is the total rotational kinetic energy of the two SIDs, 



W ; 



'r(Eg + Eg)^27rrdr 
dr 



(83) 



is the gravitational energy in the composite MSID system, 
U= I (a^E;; + a^Eg)27rrdr (84) 



is the sum of stellar and gas internal energies and 



M 



2r 



27rr dr 



(85) 



is the magnetic energy contained in the gaseous MSID com- 
ponent. 

For a full composite system with F = 1, we combine 
relations ffH) . ifTIjl . (THJ and with expressions ((H2J — JHSl 
to obtain 

T = / -Eg(rn,)^27rrdr + / -E^oi^^gf^nrdr 

Jo 2 2 ^gg^ 



at{D! + 1 



2 , 2 

+ [D^, + 1 + AV(2/3) - l//3]5 



2G(1 + (5) 



R 



expression ^ would give T/|W - M\ = 0.3281, T/jW - 
M\ = 0.1941 and T /\W - M\ = 0.0963, respectively. As 
another example, we set |m| = 2, 5 — 0.2, (3 = 2 and A'^=l. 
From equation (ITOJ, we have Dl = 0.9508, Dl = 0.1001 
and Dl = —0.4103; the third branch of solution should be 
ignored as Dl < 0. The corresponding values of T /\W — A4\ 
for the first two Dl solutions are T/|W - = 0.2700 and 
T/|W -M\= 0.0986, respectively. 

From these numerical estimates, it is clear that the ratio 
T /IW — A4\ can be much smaller than 0.14 (mainly through 
the perturbation mode of the j/2 branch), analogous to an 
unmagnetized composite SID system (Lou & Shen 2003). 
Shu et al. (2000) suggested a correspondence between T/| VV| 
ratio and the onset of secular bar-like instabilities in a sin- 
gle fluid SID. It seems natural to extend this suggestion to a 
single MSID (Lou 2002) and to an unmagnetized composite 
system of two SIDs (Lou & Shen 2003) . By this analogy, we 
now propose a further extension of this correspondence of 
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|m|=2, ^=^ 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 8. Curves of ratio T/IVV — yVf | versus S for different values 
of /3 = 1.5, 2, 2.5 for the aligned |m| = 2 case with = 1. 



aligned case 
5=0.2, (3=1.5, ?i^=0.09 



m 



Figure 9. The solution of the j/2 branch versus |m| with 
parameters S = 0.2, /3 = 1.5 and = 0.09. The smallest 
occurs at \m\ = 9. 



ratio T /\W — A4\ and the onset of secular bar- like instabil- 
ities in a composite system composed of one SID and one 
MSID. It then appears, as in the case of Lou & Shen (2003), 
the threshold of T/|VV — ratio can be considerably low- 
ered in a composite SID or MSID system. Qualitatively, this 
illustrates that the mutual gravitational coupling tends to 
make a disc system more unstable (Jog & Solomon 1984a,b; 
Bertin & Romeo 1988; Romeo 1992; Elmegreen 1995; Jog 
1996; Lou & Fan 1998b; Lou & Shen 2003; Shen & Lou 
2003). 

We now explore trends of variations for the ratio 
T/\W — M\ of the j/2 branch. Setting |m| = 2, we present 
curves of T /\W — M\ versus S when parameters f3 and A'^ 
are specified (see Fig. |SJ. As there are three solution 
branches involving three parameters 5, (3 and A^, values of 
ratio T/|VV — A4\ can be diversified. Consequently, insta- 
bility properties of a composite (M)SIDs system are more 
complex than in a single SID. Moreover, possible perturba- 
tion modes of the solution of j/2 branch can vary with 
different values of azimuthal wavenumber \m\ as shown in 
Figs. P and nn 

For parameters S = 0.2, /3 = 1.5 and A^ = 0.09, 
attains its minimum value around jm| — 9, while for param- 
eters S = 0.2, /3 = 1.2 and A^ — 3.61, the minimum value of 
Dg occurs at |m| = 2. In association with the smallest value 
of Dj, the most vulnerable stationary perturbation config- 
uration varies with chosen parameters. In other words, a 
slightly different choice of parameters may lead to entirely 
different MHD perturbation configurations. Again, this ex- 
ample shows the diversity and complexity in a composite 
MSID system and indicates that bar-like instabilities (i.e. 
m — 2) may not necessarily be the dominant instability in 
a composite MSID system. 



aligned case 
5=0.2, (3=1.2, l^=3.6^ 



m 



Figure 10. The solution of the j/2 branch versus |r?i| with 
parameters S = 0.2, /3 = 1.2 and A^ = 3.61. The smallest 
occurs at \m\ = 2 in this case. 



4 UNALIGNED MHD PERTURBATIONS OF 
LOGARITHMIC SPIRALS 

In this section, we analyze unaligned perturbation configu- 
rations of logarithmic spiral structures in a composite MSID 
system. We first derive the stationary dispersion relation for 
both full and partial composite MSID system. We then fo- 
cus on the full case {F = 1) and address the problem of 
axisymmetric marginal instabilities with \m\ = 0. 



4.1 Stationary Dispersion Relation for 
Logarithmic Spirals 

To construct stationary MHD perturbation configurations 
of unaligned logarithmic spirals in a composite MSID sys- 
tem, we pick the following potential-density pair (Syer & 
Tremaine 1996; Lynden-Bell & Lemos 1993; Shu et al. 2000; 



12 Y.-Q. Lou and Y. Zou 



Lou 2002; Lou & Fan 2002; Lou & Shen 2003), namely 

-3/2+^^ , (90) 



-3/2 + ia 



(91) 



where CTs and cjg are two constant coefficients and a is 
the radial scaling parameter related effectively to the radial 
wavenumber, together with 

(92) 



V — Vsr ' + VgV ' 



(Js and (Tg respectively by 



Vs = — 27rGA/'m(a)crs 



(93) 



Vg = -27rGAA„(a)a5 , (94) 

with A/'m (a) = K{a, m) being the Kalnajs function (Kalnajs 
1971). Consistently, we may write 

(95) 



Ug = UgT 



and 



Ua = UsV 



-1/2 + ia 



(96) 

where Ug and Us are two constant coefficients. Using expres- 
sions (|5TJ, ^ and IHSJ, equations ^ and S 
lead to 



mQgag 

^0 



+ iug ia + 



Cl 



2Q?gr'^ 2 



2j^2 2 2 2 <^2 

m S2„r — Kgr — Oa 



+ Vs + Vc 



( cl 



(97) 



- 2 



itit, = I ia + 2 — 



\2Q?gr'^ 2 



2!r22j.2 



mQgag 

^0 



(98) 



, 3 \ / o;„(T„ 
-|-mr2gr I + 2 ) 1 + Vs + Vg 



+ Vs +Vc 



i7? : 



rl/2i 



-lUgT 



Z = 



m.Q,gr 



Qgr 



iaBg 



mQgr 



(99) 



(100) 



A combination of equations 1971 and 19811 together with re- 
lations 1931 and 1941 give 



k 

c + 



^ -27rGAA„(a)j 



^0 



ag= - -^j 27rGAA™^(a)cr, 



where the three coefficients IC, A and C are defined by 
^ - ( ™ 2 ^^ ^ ^ (a^ + - l)Gi 



3/ Gi 



(101) 

(102) 
(103) 



C ; 



2rEg V2J^ir2 



(104) 



For stationary coplanar perturbations in the stellar disk in 
parallel, we set oj = in equations 1371 — 13911 and use equa- 
tions lITO — lIMIl and lISBl to obtain 

,2 



+ " + 7 ) ( ^ - 27rGM„(Q 
- — (m'-2)n?r-'' 



(105) 



Equation 11051 can also be obtained by setting Ca = in 
equation pOll and exchange subscripts s and g. 
By equations p()5|l and ll()l|l . we derive 



02 m2 



2^1 

-a +- 



X - 27rGM„(Q)E('', 



+ 



Cl + al 2-KGNm(a)T.l\ ,2,2 



2o2 



+ - 



2^2 

A 



4r2 



27rGAA™.(a)Eg\ , 2 



2n" 



: 47r^G^AC(a)EgEf, ( m'^ + + 



2 1 

" -4 



(106) 

In the absence of the gravitational coupling between the stel- 
lar SID and the gaseous MSID represented by the right-hand 
side term of equation 11061 , the left-hand side of dispersion 
relation 11061 would be reduced to two separate dispersion 
relation factors. The first dispersion relation would be 



2_L 1 

-a +- 



^-27rGAA™(Q)Ef, 



= 
"(107) 

for stationary logarithmic spirals in a single stellar SID with- 
out magnetic field. Substituting expressions 1111 and 1191 of 
Qs and Eq with F — 1 into equation 11071 . we have 



= 1 



(108) 



which is simply equation (37) of Shu et al. (2000). The sec- 
ond dispersion relation would be 



4^4 



2nl + 



27rGAA„(a)ES 



2^2 ^ rn^C'l 



27rGAA„(a)E»\ , 2 



=0 

4r2 

(109) 

for stationary logarithmic spirals in a single gaseous MSID 
with coplanar magnetic field. Equation 11091 is simply equa- 
tion (3.4.15) of Lou (2002). 
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Without magnetic field with Ca = 0, equation II1U6II 
is equivalent to the dispersion relation (86) for stationary 
logarithmic spiral configurations in a composite SID system 
of Lou & Shen (2003) as it should be. 

For numerical computations, there are two useful for- 
mulae of A/'m(a), namely, the recursion relation 

M„+i(Q)AA™(a) = [(m + 1/2)2 + 0^]"' 

and the asymptotic expansion 

AA„(a)R.(m' + Q' + l/4)-'/' 

when + > 1 (Kalnajs 1971; Shu et al. 2000). 

For the purpose of examining phase relationship be- 
tween the azimuthal magnetic field perturbation be and the 
surface mass density perturbation Eg of the gaseous MSID 
in the full case of F = 1, we derive the following relations. 
From equations 19711 and 11051 1 with F — 1, we readily obtain 



(110) 



(111) 



2TTGma„ 



[ia + Cy{2Qlr^)-3/2]{n,r) 

{1 + 5){1-Dl) AJV^ia) 



where the two coefficients A and B are defined by 
A={m^ + a'^ + 1/4) - (m^ - 2)D? 

and 



B = 



'(m' + QVl/4)(l+Df) 



(l + '5) 

It follows from expressions (irnUl and lITT^ that 
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that will be discussed later in our analysis. 



(112) 



(113) 



(114) 



(115) 



4.2 Marginal Stability for Axisymmetric 
Disturbances 

While the |m| = case is somewhat trivial as a rescaling in 
the aligned case, it is of considerable interest in the special 
'spiral' case with radial oscillations. In this subsection, we 
analyze this situation for a full composite MSID system with 
F = 1. A substitution of expressions Ulljl . dJ, H19^ and l|20|l 
into equation p()6|l with = C\/a^ leads to 
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Figure 11. The marginal stability curve of D% versus a with 
m = 0, <5 = 0.2, P = 1.5 and = 1. 
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Figure 12. The marginal stability curve of versus a with 
m = 0, <5 = 0.2, /3 = 1.5 and A^ = 3.61. 



By recursion formula 11101 and asymptotic expression lllH 
for the Kalnajs functions, we have an approximate expres- 
sion oi No{ct), namely 



AAo(q) = {a +<i/A)/{a^ + l/A)U2{,a) 



(117) 



= (a' + 9/4)/[(q' + l/4)(a' + 17/4)^/"] 

A substitution of approximate expression 111711 for Na{a) 
and expression 12311 for into dispersion relation 11161 
yields a cubic algebraic equation in terms of D^; and sim- 
ilarly, a substitution of approximate expression 11171 for 
A/'o(a) and expression 1241 for yields another cubic alge- 
braic equation in terms of . As noted earlier, we consider 
the cubic equation of because of the fact that < 
when /3 > 1. For different values of parameters 5, (3 and A^, 
we show curves of versus a in Figs. Ill( - 115l 

For the positive portions of solution DI with typical pa- 
rameters, it turns out that the basic features of a Di versus 
Q profile are qualitatively similar to those of a single SID 
[see fig. 2 of Shu et al. (2000)] and to those of a composite 
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Figure 13. The marginal stability curve of versus a with 
m = 0, <5 = 0.2, 13 = 10 and = 1. 
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Figure 15. The marginal stability curve of versus a with 
m = 0, <5 = 0.2, 0=1 and = 0.09. For (3^1, the lower ring 
fragmentation branch would gradually rise above the horizontal 
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Figure 14. The marginal stability curve of versus a with 
m = 0, 5 = 1, /3 = 10 and A^ = 1. 



SID system (see figs. 5 — 10 of Lou & Shen 2003), except 
that now trends of variation for witli tiiree parameters 
become more complicated. For example, the lower ring frag- 
mentation curve first discussed in a single MSID (Lou 2002) 
becomes negative or indiscernible in most cases (see Ap- 
pendix B). For cases of 13 around 1, however, the lower ring 
fragmentation curve of may become positive as shown in 
Fig.lTKl 

As parameters 5, f) and vary, we note several trends 
of variation in the profiles of Dl. For example, when A^ 
increases for fixed 5 and jS, the ring fragmentation curves 
of DI seem to be raised, and the collapse regime tends to 
be enlarged by comparing Figs. ITTI and 1121 In other words, 
the presence of coplanar magnetic field tends to increase 
the chance of collapse but reduce the danger of ring frag- 
mentation. This feature is qualitatively similar to the case 
of a single coplanar magnetized SID studied by Lou (2002) 
(see fig. 1 of Lou 2002). As /3 increases for fixed S and A^, 



a 



Figure 16. The marginal stability curve of versus a with 
m = 0, <5 = 0.2, /3 = 1 and A2 = 3.61. Compared with Fig. QSl 
the increase of A moves the lower ring fragmentation downward. 



the curves of Dl tend to be lowered by comparing Figs. 1111 
and 1131 When 5 increases for fixed (3 and A^, the ring frag- 
mentation curves of decrease, and the collapse regime 
shrinks by comparing Figs. 1131 and 1141 Note that as 5 in- 
creases {5 = Eq/Eq), the infiuence of gas disc on the stellar 
disc becomes more important. Therefore, the gravitational 
coupling with another disc tends to suppress the collapse 
but enhance the danger of ring fragmentation. 

One behaviour for the variation of the lower ring frag- 
mentation curve appears somewhat surprising. It seems to 
be an characteristic feature involving magnetic field (Lou 
2002). As already noted earlier, when /3 approaches 1, this 
curve could rise above the a axis to become a possible per- 
turbation mode (Fig. I15|l. However, comparing Fig. El with 
Fig. 1161 we find that in typical parameter regimes, the in- 
crease of A would force the lower ring fragmentation curve 
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to go downward. This is completely different from that of a 
single MSID (see fig. 1 of Lou 2002). 

Through numerical explorations, we realize that the lo- 
cation a of the vertical asymptote for the ring fragmentation 
branch seems to be independent of S, (3 and A^. By letting 
the coefficient of the cubic term in the cubic algebraic equa- 
tion oi y = Ds to vanish, the condition that this critical 
value Uc must satisfy is once again 

AAo(a)(a^-hl/4) =2 , (118) 

which is exactly the same as that of Shu et al. (2000), Lou 
(2002) and Lou & Shen (2003) for three related but different 
SID systems. From equation pi7^ and asymptotic expan- 
sion l|lll|l . this critical value Oc is estimated to be 1.759. 

Physical interpretations for the marginal instability 
curves are clear. As a parameter is a measure for the ra- 
dial wavenumber, a smaller a corresponds to a larger radial 
scale of perturbations. Thus, the collapse regime indicates a 
rotationally modified Jeans instability to which a compos- 
ite MSID system is vulnerable when a radial perturbation 
scale is sufficiently large (i.e. a sufficiently small a). How- 
ever, when Dg becomes larger, the conservation of angular 
momentum works against the inward gravitational force and 
the Jeans instability is suppressed. 

Nevertheless, a composite MSID system is also vulnera- 
ble to Toomre-type instability (Safronov 1960; Toomre 1964) 
associated with the ring fragmentation branch (Shu et al. 
2000). In various galactic contexts, there have been numer- 
ous studies to identify an effective Q parameter for a com- 
posite disc system (Jog & Solomon 1984a,b; Berlin & Romeo 
1988; Kennicutt 1989; Romeo 1992; Elmegreen 1995; Jog 
1996; Lou & Fan 1998b). In the case of a single SID, Shu et 
al. (2000) noted that the minimum of the ring fragmentation 
curve is closely related to the Toomre Q parameter (Toomre 
1964). In the case of a coplanar magnetized SID, Lou (2002) 
and Lou & Fan (2002) further found that the minimum of 
the upper ring fragmentation curve is tightly associated with 
the generalized MHD Qm parameter (Lou & Fan 1998a). For 
a composite system of two SIDs, Shen & Lou (2003) recently 
suggested a straightforward D— criterion to effectively deter- 
mine the axisymmetric stability against ring fragmentation 
(see Elmegreen 1995 and Jog 1996 for other proposed crite- 
ria). We therefore expect that in the present problem, the 
minimum of the upper ring fragmentation curve should be 
related a magnetic field modified D— criterion to effectively 
determine the axisymmetric stability against ring fragmen- 
tation in a composite system of one ffuid SID and one MSID. 



4.3 Stationary Logarithmic Spiral Configurations 

4-3.1 Behaviours of Di solutions 

Parallel to the study of aligned |m| > 2 cases, we substitute 
expressions ifTT^I . itTH . if^jl . and J^SJ into equation H106|l 
with A^ = C\/a^g and obtain a cubic algebraic equation of 
y = for unaligned spiral cases. One can formally derive 
three analytic solutions (see Appendix D) that would 
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Figure 17. Curves of versus a for a stationary unaligned 
logarithmic spiral with \m\ = 2, 5 = 0.2, /3 = 1.5 and = 0.09. 



be valuable in numerical MHD simulation studies of large- 
scale SID dynamics within the overall MHD density-wave 
scenario. 

To explore parameter regimes, we obtain Dg solution 
curves numerically for cases of |m| = 1, 2, 3, • • ■ and so forth. 
As expected, there are three solution curves of versus a 
given parameters |m|, 5, and /3. In certain parameter 
regimes, portions of the lower two solution curves may 
become negative and hence unphysical. We defer our dis- 
cussion of |m| = 1 case later. In Figs. I17( - I19l we show a set 
of numerical examples with \m\ =2. Through extensive nu- 
merical examples of stationary logarithmic spiral cases from 
|m| = 2 to \m\ = 5, it is clear that behaviours of Dl solu- 
tions and properties of phase relationships among perturba- 
tion variables for |m| > 2 cases are very similar to those of 
jmj = 2 logarithmic spiral. We therefore carefully examine 
the |mj = 2 spiral case for a relatively simple yet meaning- 
ful analysis. In addition, we plot curves of versus f3 by 
specifying a, 5 and A^ in Fig. 1201 

As in the aligned case, we shall use yi, j/2 and 2/3 to 
denote the upper, middle and lower solution branches, re- 
spectively. It is apparent that the three solution branches 
do not intersect with each other. Similar to the case of ax- 
isymmetric disturbances, the lowest yz solution branch does 
not appear above zero in many cases, but for cases with rel- 
atively small (3 value (e.g. /3 = 1), j/3 solution branch can be- 
come a possible stationary perturbation mode (see Figs. 1191 
and l201 . We note that always increases with increasing a 
for the upper two branches yi and 1/2 (see Figs. I17( - I19ll but 
decreases with increasing 13 (see Fig.l2()|l. The uppermost yi 
branch remains always positive, while the middle 1/2 branch 
may cross the horizontal a axis to become negative in some 
cases (see Fig. ll8ll . similar to the situation without involving 
coplanar magnetic field (Lou & Shen 2003). 

In Fig. 1201 there exist a /3ci for the middle j/2 branch and 
a (3c2 for the lowest j/3 branch for given a, 5 and A^ such that 
J/2 and 2/3 become zero and negative with increasing value 
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Figure 18. Curves of versus a for a stationary unaligned 
logarithmic spiral with \m,\ = 2, 5 = 1, /3 = 3 and = 0.09. 




Figure 20. Curves of versus /3 for a stationary unahgned 
logarithmic spiral with |m| = 2, S = 0.2, a = and = 1. 
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Figure 19. Curves of Dj versus a for a stationary unaligned 
logarithmic spiral with |m| = 2, 5 = 0.2, (3 = 1 and A^ = 1. 



of l3, respectively. Parallel to the aligned case, we examine 
behaviours of the solutions and obtain the two critical values 
of (3 analytically. For the |m| = 2 case, the analytic forms of 
Pel and /3c2 are given by 

22 



and 



Pel = ^[-B + - AACy^^] 



(119) 



(120) 



where definitions of the three coefficients A, B and C are 
given in Appendix A. When /3 < /3ci, the middle solution 
branch y2 = is positive, while for /3 > Pci, the middle 
solution branch 1/2 = becomes negative and hence un- 
physical. For the case of 5 = 0.2, A'^ = 0.09 and q = 0, we 
obtain /3ci = 2.409 from equation l|119|l . consistent with Fig. 
1171 in Fig. 1171 we have /3 = 1.5 < /3ci and the middle branch 
y2 is positive at a = 0. For A'^ = 1.00, expressions l|119|l 
and 112011 give /3ci = 2.468 and /3c2 = 1.140, respectively, 
consistent with Fig. 1201 



4-3.2 Phase relationships among perturbation variables 

For stationary logarithmic perturbation spirals, we now ex- 
amine spatial phase relationships between perturbations of 
the surface mass densities of the two SIDs (i.e. /ig and ^s), 
and between perturbations of the azimuthal magnetic field 
and the surface mass density of the gaseous MSID (i.e. bg 
and S]*). This information might provide necessary and use- 
ful clues for optical and synchrotron radio observational di- 
agnostics of spiral galaxies. 

From equation Hl()5^ . one obtains 
t._2 , _2 , T /^N^l ^^ 

Ms 



= (m^ -f + 1/4) 



r 2 , 2 
[m + a ■ 



1/4) 



jl + S) 
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- A/'m(a) 
{l + 5)y 



(121) 



1 + y 



with y = Ds- Substitution of expression 1121^ into the cubic 
equation of gives a cubic equation of fig/fj,s, and we plot 
several curves of ^ig/fis versus /? with parameters a = 1, 
S = 0.2, A^ = 0.09 and 3.61 in Figs. 1^ and 1^ 

For the phase relationship between and Z, we start 
from equation (11151 . For a sufficiently small a, correspond- 
ing to a relatively open spiral structure, one may keep linear 
terms of a and ignore the term. It follows that Z is ei- 
ther ahead of or lagging behind /ig with a phase difference of 
about ~ 7r/2. On the other hand, for a sufficiently large a in 
the tight-winding or WKBJ regime, one may drop the imag- 
inary part in equation 111511 (see Lou 2002). When accuracy 
is not that crucial in a quantitative analysis, this expression 
may also be roughly used in cases of a ~ 1. We can rewrite 
expression l|115|l approximately as 



Ms 
Z 



(1 + S){l-D] 
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where A and B are defined by equations H113^ and l|114|l . 
respectively. When parameters a, /3, S and A^ are specified 
in equation 112211 . solutions of y = correspond to val- 
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Figure 21. Curves of fig/fJ.s versus fS for a stationary logarithmic 
spiral with |m| = 2, S = 0.2, a = l and = 0.09. 
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Figure 23. Curves of real parts of fig/Z versus /3 for a stationary 
logarithmic spiral of |m| = 2, S = 0.2, a = 1 and = 0.09. 
Imaginary parts of Hg/Z are not shown here. While one branch 
of real part of fJ.g/Z corresponding to yi remains always negative, 
this branch does not remain always lower than the real part of 
Hg/Z branch corresponding to j/2. 
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Figure 22. Curves of /J^g/fJ-s versus fS for a stationary logarithmic 
spiral |m| = 2, 5 = 0.2, a = 1 and A-^ = 3.61. The uppermost 
branch corresponding to j/3 = is only discernible for /3 < 1 
and is not shown here. The branch corresponding to 1/2 = 
may cross the horizontal /3 axis. 



ues of i-ig/Z. We then present curves of j-ig/Z versus /3 with 
parameters a = 1, 3 = 0.2, — 0.09 and 3.61 in Figs. 
1231 and 1^ By equation I12H . we note again that a smaller 
Hg/fJ-s corresponds to a larger y = D^, reminiscent of the 
same variation trend of correspondence in the aligned case 
described in Section 3. Such a relatively simple correspon- 
dence does not hold for fJ,g/Z that should be analyzed with 
extra care. Through extensive numerical explorations, we 
note that the highest branch of Hg/Z always relates to the 
lowest j/3 branch, while the j/i branch of Hg/Z (being always 
negative) does not remain always less than the 2/2 branch of 
IJ.g/Z as shown in Figs. l23l and l24l These phase relationships 
among branches of j-ig/fis and ^ig/Z and the corresponding 
branches of y = are shown explicitly in Figs. I21H24I 

The unaligned spiral cases here are fairly similar to the 
aligned cases of Section 3, except that the imaginary part of 
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Figure 24. Curves of real parts of fJ,g/Z versus /9 for a stationary 
logarithmic spiral of \m\ = 2, 5 = 0.2, a = I and A^ = 3.61. 
Imaginary parts of fJ.g/Z are not shown here. The uppermost 
branch for real part of Hg/Z is only discernible for /3 < 1 and is 
not shown here. The branch of real part of fJ.g/Z corresponding 
to 3/2 = may cross the horizontal /3 axis. 



I-Lg/Z does not vanish in general. For the stationary mode of 
yi branch, the phase relationship between the surface mass 
density of the two SIDs (i.e. fig/fis) and the phase relation- 
ship between the surface mass density of the MSID and the 
azimuthal magnetic field (i.e. i-ig/Z) are always out of phase 
(see Figs.l^Tl-l^.t 



I To be more precise, fJ,g/Z being 'out of phase' means that the 
real part of /jg /Z is negative. As a ranges from sufficiently small 
to relatively large values, the phase difference between fig and Z 
ranges from 7r/2 to tt, either ahead or lagging behind. Likewise, 
lig/Z being 'in phase' implies that the phase difference between 
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For the unaligned perturbation mode of j/2 branch in 
most cases, ^ig/fis is in phase, while the real part of fig/Z is 
always negative. Taking into account of the imaginary part 
of Hg/Z, the phase difference between fig and Z is ir /2 ^ n 
(see Figs. l21l and l23^ . This j/2 branch has a zero point at /3co 
(where /ig — 0) that increases with increasing A^. When 
becomes sufficiently large (e.g. A^ = 3.6f) with (3 > 1, the 
value of /3cO can become larger than f. ft is then possible 
for fig/ fj,s and fig/Z to become zero at /3 = /3co and reverse 
their respective signs for 1 < /3 < /3cO as shown in Figs. 1221 
and 1241 Taking into account of the imaginary part of fig/Z, 
this implies that the phase difference of fig and Z switch 
from between 7r/2 ~ tt to between ~ 7r/2 as /3 varies from 
f} > Pco to /3 < Pco- Meanwhile during the same (3 transition, 
the phase difference between fis and fig switch from in phase 
to out of phase. The special situation of /Xg = at /3 = /3cO, 
occurring for stationary perturbation mode of y2 branch, is 
very similar to that of the aligned case studied earlier. With 
fig = Jg =G but fis i= 0, Z ^ Q, ^g ^ Q smA Ug / 0, thls 
may correspond to a unique solution. 

For the stationary perturbation mode of ya branch, both 
fig/ fis and the real part oi figjZ remain positive and there- 
fore fig and fis are always in phase (see Figs. l2f H24|l : taking 
into account of the imaginary part of fig/Z, the phase dif- 
ference between fig and Z ranges from to 7r/2. 

The analytic expression of Pco for |m| = 2 and a = 1 is 
given by 

1 



609(1 + 5) ' 



[A^(17V21 4-685 -h( 



(123) 



-fl68 5 +42 V21 + 168 



that increases with increasing A^. For 5 — 0.2 and A^ = 0.09 
in equation 1123II . one has /3co = 0.5589 that is not shown 
explicitly in Figs. 1211 and 1231 (see Fig. El in Appendix E). 
For 5 — 0.2 and A^ = 3.61, one has /3co — 1.3271, consistent 
with Figs. 123 and El 



4-3.3 Solution behaviours of unaligned \m\ = 1 case 

We now turn to the somewhat special case of |m| = 1. 
Mathematically, the qualitative difference between cases of 
|mj = 1 and |m| > 2 is perhaps due to the fact that 
Nmiot) < 1 for |m| > 2 but otherwise for \m\ — 1 (Lou 
& Shen 2003). We now examine behaviours oi y = of the 
|mj = 1 case shown in Figs. I25l and 1261 

ft turns out that solution curves of |m| = 1 case indeed 
differ from those of |m| = 2 case. For example, the middle 
and lower branches of y = Dg remain always negative and 
hence unphysical (see Figs. I25jandl26^. Portions of the up- 
per branch of j/i = may become zero or even negative 
for a sufBciently large /3 (Fig. 1251 . That is, when /3 exceeds 
a critical value /3c3, yi = would become negative. The 
analytical form of this critical /3c3 is given expression 112411 . 



fig and Z ranges from to 7r/2 as a varies from relatively large 
to sufficiently small values. 
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Figure 25. Curves of y = versus /3 for stationary unaligned 
perturbations of |m| = 1 with 5 = 1, a = 1 and = 3.61. 
Solution branch of yi = may go across the horizontal /3 axis 
and has a root at about 56. 
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Figure 26. Curves of j/ = versus a for stationary unaligned 
perturbations of \m,\ = 1 with S = 0.2, /3 = 1.5 and A^ = 1. 



These properties are similar to those in the case of a com- 
posite unmagnetized SID system studied by Lou & Shen 
(2003). The three solution branches of t/ = do not inter- 
sect with each other. The increase of /3 moves both yi and 2/2 
branches downward (see Fig. 1251 and slightly lifts the lowest 
j/3 branch. 

One can derive an expression for the critical point of (3 
where yi = 0. For a = 1, the analytical form of this /3c3 is 

= [1 + 21s + 5A^ + 20A^5 



" (125-8)' 
+(81 + 815^ + 430A^5 + 1625 + 300A^5^ 
+130A^ + 140A*5 + 580A*5^ - 15A*)'''^1 



(124) 



For 5 = 1 and A'^ — 3.61, expression 1124^ gives f3c3 ~ 
56.143, consistent with the upper curve in Fig. 1251 There is 
no essential mathematical difficulty of obtaining a more gen- 
eral analytical form of Pc3 for arbitrary values of a, although 
the expression may appear somewhat complicated. 
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Figure 27. Solutions of for tlie 1/2 brancti versus \m\ with 
parameters a = 0, S = 0.2, (3 = 1.5 and = 0.09. The smallest 
appears at \m\ = 10. 
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Figure 28. Solutions of for the y2 branch versus |m| with 
parameters a = 0, 6 = 0.2, /? = 1.2 and = 3.61. The smallest 
appears at |m| = 2. 



Physically, one may exchange the role of |m| and a as 
they both contribute to the total dimensionless wavenum- 
ber. Given a fixed a value, a composite MSID system may 
support different values of |m| (see fig. 3 of Shu et al. 2000 
and figs. 14 and 15 of Lou & Shen 2003). This is shown 
here in Figs. I27| and 1281 Analogous to the aligned case, the 
smallest value of j/2 = occurs at an azimuthal periodic- 
ity of |m| = 10 with parameters a = 0, S = 0.2, l3 = 1.5 
and — 0.09. As a comparison, for the same S — 0.2 and 
a = yet a smaller P — 1.2 and a stronger magnetic field 
= 3.61, the smallest value of 1/2 = Dl then occurs at 
|m| = 2. Again, we see that the stability of a composite 
(M)SID system appears to be more complicated than those 
of a single (M)SID system by the involvement of three pa- 
rameters. As parameters 5, f3 and A^ vary, vulnerable con- 
figurations with the smallest value change accordingly. 



5 DISCUSSION AND SUMMARY 

In this final section we outline some potential applications 
of our analysis in the context of a nearby spiral galaxy NGC 
6946, and then summarize our theoretical results and sug- 
gestions. 



5.1 Galactic Applications 

By the hydrodynamic density-wave scenario for a disc 
galaxy, large-scale high gas density spiral arms are vulner- 
able to molecular cloud and star formation activities on 
smaller scales and are therefore luminous in optical bands. 
As synchrotron radio observations also reveal large-scale spi- 
ral structures in association with optical spiral arms (e.g. 
Mathewson et al. 1972; Sofue et al. 1986; Beck et al. 1996), 
it is necessary to incorporate galactic magnetic field em- 
bedded in the interstellar medium (ISM) and an relativistic 
cosmic-ray gas (Lou & Fan 2003) in the MHD density-wave 
scenario for a deeper physical understanding (Fan & Lou 
1996; Lou & Fan 1998a). 

The interlacing of optical and magnetic spiral arms in 
the nearby spiral galaxy NGC 6946 was first revealed by 
Beck & Hoernes (1996) in the inner disc portion with an al- 
most rigid rotation. Similar features of displaced optical and 
radio arms were also reported in portions of spiral galax- 
ies IC 342 (e.g. Krause et al. 1989) and M83 (NGC 5236; 
e.g. Sukumar & Allen 1989; Neininger 1992; Neininger et al. 
1993). These observational results led Fan & Lou (1996) to 
introduce theoretical concepts of slow MHD density waves 
(SMDWs) and fast MHD density waves (FMDWs) in mag- 
netized spiral galaxies (Lou & Fan 1998a; Lou, Yuan & Fan 
2001) and to show that perturbation enhancements of az- 
imuthal magnetic field and gas mass density of SMDW are 
phase shifted by > n/2 relative to each other. However, 
several wavelet analysis on multi-wavelength data of NGC 
6946 (e.g. Frick et al. 2000, 2001) further revealed large-scale 
spiral arms extended well into the outer disc portion with a 
more or less fiat rotation curve (e.g. Tacconi & Young 1989; 
Sofue 1996; Ferguson et al. 1998), that seems to challenge 
earlier arguments that SMDWs are largely confined within 
the inner disc portion of almost rigid rotation. 

For coplanar MHD perturbation structures of station- 
ary logarithmic spirals in an MSID model with a flat rota- 
tion curve, Lou (2002) has shown the existence of SMDWs. 
For stationary SMDWs of Lou (2002), logarithmic spiral en- 
hancements of gas surface mass density and azimuthal mag- 
netic held are phase shifted by about 7r/2 ~ vr. As effects of 
long-range self-gravity, magnetic field and disc differential 
rotation are included in the model analysis, SMDW patterns 
in a single MSID model can indeed support extended man- 
ifestations of interlaced optical and magnetic spiral arms 
that persist well into the outer disc portion with a flat rota- 
tion curve (Lou & Fan 2002). Moreover, Lou (2002) further 
suggested that significant phase shifts between ~ 7r/2 re- 
lated to stationary logarithmic spiral structures of FMDWs 
within a single MSID. 
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Given idealizations of our composite MSID model, 
the three possible stationary logarithmic spiral patterns 
might be conceptually relevant to magnetized spiral galax- 
ies. Specifically, for the upper two solution branches of yi 
and y2, the phase difference of the stationary logarithmic 
spiral enhancements between azimuthal magnetic field and 
surface mass gas density is between 7r/2 ~ tt, supporting the 
idea that interlaced optical and magnetic spiral arms may 
persist through the inner disc portion to the outer portion 
(with a phase difference > 7r/2). For a sufficiently large a 
in the tight-winding or WKBJ regime, this phase difference 
approaches ~ n, while for a sufficiently small a in the open 
regime, this phase difference approaches 7r/2. Note that in 
our modelling, the phase difference of ^ n/2 occurs mainly 
through the 3/3 solution branch that is usually unphysical 
with being negative. 

As sketched earlier (e.g. Lou & Fan f998a, 2000a, b; 
Lou et al. 2002), a typical disc galaxy involves a stellar disc, 
a magnetized gas disc and a massive dark matter halo. Our 
theoretical model of a composite MSID system, idealized 
and simplified in many aspects, does contain these basic el- 
ements. With necessary qualifications in mind, it is then of 
considerable interest to reveal possible types of large-scale 
structures and their characteristic features in the station- 
ary MIfD density-wave scenario and relate them to galactic 
diagnostics and observations. 

For the t/i solution with largest value of D^, station- 
ary logarithmic spiral density arms in gaseous MSID and in 
stellar SID are spatially out of phase; in particular, spiral 
gas density and magnetic field arms are siginificantly phase 
shifted relative to each other with a phase difference between 
■k/2 ~ n. 

For the 7/2 solution with middle value of D^, stationary 
logarithmic spiral density arms in gaseous MSID and in stel- 
lar SID are spatially in phase; meanwhile, spiral gas density 
and magnetic field arms are siginificantly phase shifted rel- 
ative to each other with a phase difference between 7r/2 ~ tt 
for /3 > f and (3 > /3co. At the qualitative level, this case 
appears comparable to interlaced optical and synchrotron 
radio arms of NGC 6946 (Beck & Hoernes 1996; Fan & Lou 
1996; Frick et al. 2000, 2001; Lou & Fan 2002). On the other 
hand, for 1 < (3 < jSco, stationary logarithmic spiral density 
arms in gaseous MSID and in stellar SID are spatially out 
of phase; meanwhile, spiral gas density and magnetic field 
arms are phase shifted relative to each other with a phase 
difference between ~ 7r/2. 

By the above analogy of unaligned spiral configurations, 
our results for aligned stationary MHD perturbations with 
\m\ = 2 may bear possible diagnostics for optical and syn- 
chrotron radio observations of galactic bar structures and 
their spatial phase relationships. It is straightforward to 
compare Figs. |11-J7| of the aligned cases with Figs. I21H24I 
and infer spatial phase relationships of stellar bars, gas bars 
and synchroton radio bars involving magnetic fields. For ex- 
ample, by Figs. 0] and 1^ the yi solution describes a config- 
uration in which the stellar density bar and the gas density 
bar are phase shifted relative to each other with a n phase 



difference, while the gas density bar and the synchrotron ra- 
dio bar (i.e. magnetic bar) are phase shifted relative to each 
other with a tt phase difference. The 2/2 solution describes 
a configuration in which the stellar density bar and the gas 
density bar coincide, while the gas density bar and the syn- 
chrotron radio bar (i.e. magnetic bar) are phase shifted rel- 
ative to each other with a tt phase difference. By Figs. |K| 
and|7| the y\ solution describes a qualitatively similar bar 
configuration as the j/i configuration of Figs. 2] and |H] For 
P > Pco [see definition 17711 for /3co], the J/2 solution describes 
a qualitatively similar bar configuration as the 1/2 configu- 
ration of Figs. 2] and |S| For 1 < /3 < /3cO, the 2/2 solution 
describes a configuration in which the stellar density bar 
and the gas density bar are phase shifted relative to each 
other with a tt phase difference, while the gas density bar 
and the synchrotron radio bar (i.e. magnetic bar) coincide. 
In context of observations, the prominent bar structure of 
the nearby barred spiral galaxy M83 (e.g. Sukumar & Allen 
1989; Neininger 1992; Neininger et al. 1993) is worth being 
analyzed in more details in both optical and synchrotron 
radio bands. 

5.2 Summary 

In terms of theoretical development of MSID model and in 
reference to previous theoretical analyses of Shu et al. (2000) 
on zero-frequency (i.e. stationary) aligned and unaligned 
perturbation configurations of an isopedically magnetized 
SID, of Lou (2002) on stationary aligned and unaligned per- 
turbation configurations of a coplanarly magnetized SID, 
and of Lou & Shen (2003) on stationary aligned and un- 
aligned perturbation structures in a composite system of 
two-fluid SIDs, we have constructed analytically in this pa- 
per stationary aligned and unaligned perturbation configu- 
rations in either full or partial composite system consisting 
of a stellar SID and a gaseous MSID. While this composite 
model of one SID and one MSID is highly idealized in many 
aspects, it does contain several necessary and more realistic 
elements that are pertinent to structures and dynamics of 
disc galaxies. Galactic structures may not be stationary in 
general, yet the stationary global MHD perturbation struc- 
tures might represent a subclass, provide insights for a disc 
galaxy, and serve as benchmarks for numerical MHD simu- 
lations. Given our model specifications, we have reached the 
following conclusions and suggestions. 

For aligned MHD perturbation structures, we derived 
stationary dispersion relation 1621 for both full and par- 
tial composite system of one SID and one MSID. Mathe- 
matically, there are three possible values for the square 
of dimensionless rotation parameter Ds of the stellar SID 
(or equivalently, -Dg). Physically, only those solutions with 
> and Dg > may be conceptually applicable. 

In the aligned case of barred configurations with \ra\ = 2 
and F = 1 for a full composite MSID system, the upper- 
most branch remains always physically possible, with spa- 
tial phase relationships jJ.g/ fis < and Hg/Z < (where Z 
is the amplitude of azimuthal magnetic field perturbation); 
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that is, both fj,g and fia pair and fig and Z pair are out 
of phase. These qualitative features for possible bars may 
be detectable observationally. When parameter /3 is smaller 
than a critical value Pci defined by l|71|l . the middle solution 
branch of y = is positive and thus physically possible. 
Corresponding phase relationships of fig/ fia and of fig/Z are 
in phase and out of phase, respectively. But when /3 < /3cO 
defined by expression 17711 . fig/ fis and of fig/Z can become 
zero or even change to opposite sign (see Figs. |S] and |7J. This 
special fig = case might be an additional possibility. The 
lowest branch of y = is physically possible only when 
P < Pc2 defined by expression (1721 . In many cases /3c2 is 
smaller than 1, this branch is thus invalid. However in some 
cases (see Fig.|5J, this branch can rise up to become positive 
and thus physically possible. The corresponding phase rela- 
tionships of fig/ fia and fig/Z are both always positive and 
thus in phase. 

For aligned perturbations, the case of |m| = can be 
made to effectively rescale the axisymmetric background. 
Eccentric |m| = 1 displacements may occur for arbitrary 
Ds values in a full composite MSID system. In contrast, 
for a partial composite MSID system, Dl can no longer be 
arbitrary for jml = 1 (Lou 2002; Lou & Shen 2003). 

For coplanar MHD perturbation structures of unaligned 
logarithmic spirals, we obtained the stationary dispersion 
relation H106|l for both full and partial composite MSID sys- 
tem. Mathematically, there are three possible sets of solu- 
tions for or equivalently D^. 

For stationary logarithmic spiral configurations with 
|m| = 2 and _F = 1 in a full composite MSID system, the 
uppermost Dl solution branch is always physically possible, 
with corresponding fig/ fis being always out of phase and 
with the phase difference between fig and Z in the range 
7r/2 ~ TT. The middle branch of solution is physically 
possible when (5 < jSci defined by expression 111911 . There 
also exists a special Pco defined by expression 11231 : when 
/3 decreases below f3co, the corresponding fig/fis changes 
from being in phase to being out of phase (see Fig. 1221 . 
and the phase difference between fig and Z changes from 
being within 7r/2 ~ tt to being within ~ 7r/2 (see Fig. 
1241 . The special fig — case may occur as well. The low- 
est branch of solution can be physically possible when 
/3 < /3c2 defined by expression l|120|l . Similar to the aligned 
case, while this lowest solution branch is often invalid 
as /3c2 is usually less than 1, it may become physically pos- 
sible for (3 within interval [l,/3c2] (see Figs. \W\ and 1201 for 
cases with /3c2 > 1. Correspondingly, fig/fis remains always 
in phase and the phase difference between fig and Z falls in 
the range ~ 7r/2. 

For the marginal case of |m| = with radial propa- 
gations, the collapse and ring fragmentation regimes exist, 
while the lower ring fragmentation curve, which is a novel 
feature in a single MSID (Lou 2002), does not often ap- 
pear to be positive. For cases with /3 ~ 1, this lower ring 
fragmentation branch may rise to become positive and will 
descend with increasing as shown in Figs. 1151 and 1161 
This qualitative trend differs from that of a single MSID 



studied by Lou (2002) . We further examined trends of vari- 
ation as parameters 5, jS and vary, and concluded that 
the mutual gravitational coupling between the two SIDs ap- 
pears to reduce collapse regime but to increase the danger of 
ring fragmentation. Physical interpretations on collapse and 
ring fragmentation regimes are provided (Shu et al. 2000; 
Lou 2002; Lou & Shen 2003; Shen & Lou 2003). Consistent 
with prior analyses, we realize that the ring fragmentation 
curve has a fixed vertical asymptote at A/'o(a)(a'^ + 1/4) =2 
with a ~ 1.759 approximately; this condition for the verti- 
cal asymptote remains exactly the same as that of Shu et 
al. (2000), Lou (2002) and Lou & Shen (2003). In reference 
to these previous analyses, it is highly suggestive that the 
effective Q parameter in a composite MSID system should 
be closely related to the minimum of the ring fragmentation 
curve (Shen & Lou 2003). 

The unaligned spiral case of |m| — 1 is somewhat spe- 
cial. Being negative, the two lower solution branches of 
are both unphysical. The uppermost solution branch is 
physically possible when /3 < /3c3. The analytic expression 
of /3c3 with a = 1 is given by equation 11241 . 

We have also derived the MHD virial theorem in a com- 
posite MSID system, following the suggestion of Lou (2002) 
that ratio T /|W — M\ may play the role of ratio T /|W| in 
determining the onset criterion for bar-like instability when 
coplanar magnetic field is involved. Similar to a composite 
system of two fiuid SIDs of Lou & Shen (2003), we found 
that solution of the middle y2 branch may correspond 
to T /\W — A4 \ ratios being considerably lower than the oft- 
quoted value of ~ 0.14 (Ostriker & Peebles 1973; Binney & 
Tremaine 1987)). This indicates that a composite (M)SID 
system tends to be less stable. This is just one perspective 
of showing that the stability in a composite MSID system is 
far more complicated than that of a single MSID. In partic- 
ular, additional solution branches allowed by more dynamic 
freedoms as well as more parameters 5, (3 and A'^ do intro- 
duce extra dimensions for the stability problems. For exam- 
ple, as these dimensionless parameters vary, most vulnerable 
configurations related to the smallest change. 

While our model for a composite system of MSIDs is 
highly idealized, theoretical results obtained here provide a 
conceptual basis and useful clues for optical and synchrotron 
radio observational diagnostics of spiral galaxies of stars and 
gaseous ISM with magnetic field coupled by mutual gravity 
on large scales. In particular, our analytical solutions for 
stationary MHD perturbation configurations in a composite 
MSID system as well as analytical expressions for critical 
points are important and valuable for benchmarking numer- 
ical codes designed for large-scale dynamics of magnetized 
gaseous discs. 
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APPENDIX A: 

For the spiral case of |mj = 2, a substitution of expressions 
llTTl l. lO, lO, l|2SJ and ^ into equation lUTO would lead 
to a cubic algebraic equation in terms oi y = D^. By setting 
the constant term of this cubic equation equal to zero, we 
derive conditions for y — 0, namely 

Af3'^ + Bf3 + C ^0 , (Al) 

where the three coefficients A, B and C are defined by 

^ = -578 AAa (a) 5 - 32 a'^A/'s; (a) 5 - 272 Ne {a) 5 a 
+272 N'2 (a) -I- 64 aNz (a) - 272 - 64 a 

-MS a"^ - 2725 , (A2) 

B = -400 aN2 {a) + 272 Neia) 5 a + 32 a* 

+32 a^N2 (a) (5 + 32 a" + 578 Ne (a) 5 - 34 A^a (a) 
-1122A/'2(a) +400a^ + 1122 + 144 5 A^a^ + 32 A^ 
-f34 5A^ - 144A/'2(a) - 32a''7Va(a) + 1122 J 

-1-16 aN2 {a) + 400 5a +?,2a^5 + 144 
+ 120N2{cy)5X'a +22lN2{cy)5X' -■i2a^N2{a) A^ 

(A3) 

and 

C = -850 - 264 - 264 5 - 32 a"* -f 96 A*5 a 
+96 \*a + 16 a^\'^5 - 48 a^5 \^ - 850 5 - 32 a*5 
+16 A^q" - 336 a - 255 A^ - 255 5 A^ 
-336 5 a +119 X'^5 + 119 A" - 48 a" A^ 
+850AAa(a) + 264A/'2(q) A^a^ - 96 A^AAj; (a) a^ 
+255AA2(q) A^ - 119A''AA2(q) +336a^7Vi;(a) 
+32 a'^N2 (q) - 16 a^\'^N2 {a) + 48 a*N2 {a) . 

(A4) 
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Figure Bl. The marginal stability curve of versus a with 
parameters |m| = 0, <5 = 0.2, (3 = 1.5 and = 1. Fig. IHU 
presents corresponding curves of versus a. 
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Figure B2. The marginal stability curve of Dg versus a with 
parameters \m\ = 0, S = 0.2, /3 = 1.5, = 1. Fig. IBll presents 
corresponding curves of Dg versus a. 

APPENDIX B: 

In the presence of gravitational coupling between one SID 
and one MSID, we expect three sets of independent mathe- 
matical solutions in general. However, for cases of marginal 
stabilities for spiral configurations with |m| = 0, the solu- 
tion curves qualitatively resemble those determined for a 
single SID (Shu et al. 2000) in most cases, and the two ex- 
tra solutions turn out to be negative. Moreover, the lower 
ring fragmentation curve in a single MSID (see fig. 1 of Lou 
2002) remains often indiscernible (see Figs. I11H15II . Here, 
we provide a complete structure for solution curves in Figs. 
lBTlandlB2l 

We note that the basic positive Dg versus a profile is 
qualitatively similar to fig. 1 of Lou (2002) (see Fig. IBlfl . 
However, the curve corresponding to the lower ring frag- 
mentation curve of Dg is negative in most cases, indicating 



that this branch is physically invalid (see Figs. IBl1 and lB2ll . 
On the other hand, when /3 approaches 1, this branch of D^ 
can be raised above zero and thus becomes physically valid 
in this limit (see Fig. 1151 . 

APPENDIX C: 

In this Appendix C, we provide the complete form of the 
cubic equation of y = D^ for the aligned case of coplanar 
MHD perturations in a composite MSID system for inter- 
ested readers or potential users in numerical MHD simula- 
tions. From dispersion relation (1621 and relation 12311 . we 
readily derive 

Ay"^ + By^ +Cy + V = , (CI) 
where the four coefficients A, B, C and D are defined by 

A = -16 lm|/3^ - 4m*f3'^{l + 5) + 16 \m\ 0^ 

-4,m'^\m\l3'^ -l&\m\50^ + Srn5l3^ (C2) 
-4 \m\5p'^ + 9. rn 0^ + 1& in \m\ 5 , 

B = -Am^Sp^ - l&TU P - Atu^ \m\ 0^ + 2ni^5 \^ P 
-32 \m.\ 0^ 32 |m| 5/3 - 4m* |mj 5 p"^ 
+32 \m\P-AQm^\ni\5p + A ra& /3 -f 16 rni 0^ 

+24m^P^ + 24m^\m\p'^ -32\m\5p^ (C3) 
-f 12 m'' |m| (5 /3 - 8 m'^S P - 12 m^0^ 
+45 P + 12 m'^P + 24 \m\ 5 P'^ 
-&m^5 \^ P + 12m'^ \m\ P - 4Qfn \m\ P , 

C = 24m^ \m\5 + 2\^rn -S.m'^ + 24rn \m\ + m'^A"' 
— A*m^ — 16 |m| 5 + ni^X'^ \m\ — 32 rn^ P + 24 0^ 
—8 m'* \rn\ — 16 \rn\ + 6 m^X^ \m\ 5—2 to^}? \m\ 5 
— 3 |m| A'*m^<5 + 6 m^A^ |m| — 3m^A''|m| 
— 2m''A^ jmj — 2m* A^ — 4 |m| A^ — 4 |m| 5 X^ 
-4 ■rn5 X^P - 8 m^5 P + 4 m*5 + 8 m^5 
-8 m* Imj (5 -f 8 m^ - 24 m^ jmj /3 -f 24 m*/3 
-12m*/3^ + 4m*jmj/3^ -24m^jmj5/3 + m*A*jm|5 
+4m*jmj<5/3^ +4(5A^/3- 16 jmj + 32 \m\ P 
+2 A* jmj - 16 \m\5p'^ + 32 jmj <5 /3 + 2 A* jmj 5 

(C4) 

and 

T> = —8 m^ jmj <5 -I- 2 A^m^ — 8 m* — 8 m^ jmj + m^A* 
— A*m^ — m*A* jmj — 16 P + 8 rn^ p'^ + 8 m* jmj 
— 2 m^A^ jmj 5 + 2 m*A^ jmj <5 + jmj A*m^<5 — 2 m^A^ jmj 
-f-m^ A* jmj + 2 m* A^ jmj - 2 m* A^ - 4 m*<5 P + 2 m^5 X^P 
-2 m*(5 X^P + 4 m*5 + 8 m* jmj 5 + 8 m^ - 8 m^ jmj 

-Hl6m^jmj/3 + 12m*/3- 12 m* jmj /3 - 4 m*/3^ 
-H4m* jmj P'^ - 12m* jmj (5/3+ 16m^ jmj 5/3 
— m*A* jmj 5 — 8 rn^ jmj 5/3^ +4 m* |m| 5 /3^ . 

(C5) 
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One can solve cubic equation (CI) for possible aligned MHD 
perturbation configurations in a composite MSID system by 
specifying parameters m, 5, 13 and A. 



APPENDIX D: 

In this Appendix D, we provide the complete form of the 
cubic equation of y = for the unaligned spiral case of 
coplanar MHD perturbations in a composite MSID system 
for interested readers or potential users in numerical MHD 
simulations. From dispersion relation 1621 and relation I23II . 
we readily derive 

Ay"" +By^ +Cy + V = , (Dl) 

where the four coefficients A, B, C and T> are defined by 

A = i2m^Nm{a)a^P'^ +2,2n?Nm{ct)5a^P'^ 
-l&Nm{a)l3'^ + 129.5 - I28m^/3V32m'*/3^ + 128/3^ 

-I67V„(a)(5/3^ - 64AA„(a)a^/3^ - 56AA,„(a)m^/3^ 
+32m^AA™(a)/3^ - 64AA„(a)(5Q^/3^ + 32m^AA„(a)(5/3^ 
-^&Nm{ct)5rn +?,2 m'^ 5 0^ -12?, rn 5 , 

(D2) 



C = 112 - 96 5 A^a^ - AA,„(q)A^ + A'*7V™(a) 
+16 m*A^ - 36 A^ + 20 X'^m'^S 
+48 A^a^ + 8 5 A^ - 8 m*\'^S 
+ 16m*'(5 A^ + 20m^A^ + 16 m*'A/'m(a)A^ 
+8 AC(a)A^a^ - 64a^ - 96 A^a^ 
-40m^AC(a) - 64a^(5 - 48aW™(a) 
+32 Q^A"* + 2 A*'A/;„(a)m^ - 8 A^^ 
+ 16 m'^jVm{a)5 X^a"^ P - 16m'^ X'^S - 24m^A''AA„(a)Q^ 
+48 m'^S X^a + 32 rna + 32 rn^aS 
-8m*A* + 64m^A/'m(a)A^Q^ - &-irnNrr,{a)X^a 13 
+16a*AA™(Q)(5 A^/3 + 60AA™,(q)/3 - 256m^7V„(Q)a^/3 
+80 A/'™(a)5 a 13 + 32 a''AC(a)5 /3 + 32 m^AA„(a)<5 

+96 m^a^5 /3 + 72 Mm {q)5 rn 13 - 48 a 13 
-270 5 /3 + 192 Ar^(Q)m^/3 + 224 A/'m(a)a^/3 
-48 a^(5 /3 + 96 - 192 m^A/"™ (a)/3 

+18A/'m(a)5/3 + 32a^A''(5 + 18 A^/3 



Z3 = -16m''A/'„(a)(5A^/3- 16m^7V™,(a)5A^a^/3 
-32m^AA„(a)A^a^/3 
+307V;„(a)/3 - 128m^7V™(a)Q^/3 + 64AA™(a)(5 
-32m''A/'m(a)(5/3 - 32 rn Nm{a)5 a 13 - 'i2 rna5l3 
+56 A/'m(a)5 + 64 - 240 5 /3 
-96m''/3 + 96A/™(Q)m^/3 + 112A/;i(Q)a^/3 
+64 a 5 /3 - 32 - 96 ■m^Nm{a)l3 

+l&Nm(a)5l3 + 16 A^/3 - 32m^5A^a^/3 
+64 5 X^aP - 24 A/'m(Q)(5 A^/3 - 32 m^A^a^/3 
-8 A/'™(a)A^m^/3 - lQNm{a)X^o?p - 32 a*AA„(Q)A^/3 
-2 A/'m(a)A^/3 + 64 A^a^/3 - 8 X^m^p 
-32 a''AA„(a)/? + 32 A/'„(a)5 A^ + 44 AA„(a)(5 A^ 
-144A/'m(a)(5a^/3^ - 128A/'m(a)<5m^/3^ - ?,2rna5p'^ 
+<d6m^Nm(a)aP'^ + 32 m''/?^ - 48A/'™(a)/3^ + 64a^/3^ 
+272(5/3^ - 168AA„(a)m^/3^ - 192 AA„(a)a^/3^ 
+Ua5p'^ - 32m^a^/3^ +96m''7V,„(a)/3^ 
-34AA„(q)5/3^ + 32m'*AA„(Q)(5/3^ - 32Q''7V;„(a)(5/3^ 
-200 - 240 /3 + 272 

-96 m*(5 /3 - 8 m^(5 A^/3 + 312 /3 
+32 m^S p"^ - 200 m^S P^ + 312 m^P + 16SX^P , 

(D3) 



+80 5A^a2/3- 27AA„(a)5A^/3 + 32 q-^^A^/? 
-16A/'m(a)A^m^/3 - 32 AA„(a)A^a^/3 - 64a''AA„(a)A^/3 
-4 AA„(a)A^/3 + 80 A^a^/3 + 32 a'^X'^P 
+ 112 5 + 32a''/3 + 32 a" (5/3 
-64 a'^Af,^{a)P + 24 A/"™ (a)5 A^ - 96 m^AC (a)<5 
+367V,„(a)5m^/3A^ - 96AA™,(q)(5 - 88 AA„(a)5 

-64 m^a^(5 + 96 m^TV™ (0)0^/?^ - 32 m**/?^ 
-48AA„(a)/3^ + 128 + 160 Sp^ 
-168A/'m(a)m^/3^ - 192 A/'m(Q)a^/3^ + 128a'^dp^ 
-64 m^a^P'^ +96 m'^Afrr,{a)p'^ -20 Afrr,{a)dp'^ 
-1S4 m^S - 32 m'^J\f,n{a)S p'^ - 64a'^Afm{a)5 P^ 
-16rnP^ - 184m^ +8A^ - 270 P 
+4&a^Nm{a)X^ + 160 P^ - Sm^X^Nmia) - 36A^m^ 

+64 m'' - 8 A'' +64 m''^ 
— 16 a*A'*A/'m(a) + 06rin?Hm{a)cP — 16 X*m?a^ 
+216 rnS /3 - 32 m^S P^ -16rn5 P^ 
+216 rnP + 2,2 a^Nm (a) + 64 m^Nm {a) 
+ l'g,&X^P-UNm{a) , 
(D4) 
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— 16 m*A^ — 2 A*m^5 — 64 m^A^o?^ 
+(5A^ + 8m^A^(5-16m''5A^ 
-2 m^A* - 48 a^S A^ + 16 m*Um(a)X^ 
+ 16m''Mn(a)<5 A^/3 + 87V„(a)A^a^ + 48 
-8 A^a^ - 40 rnNra (a) + 48 
-48a^A/'m(a) + 2 A''AA„(a)m^ - A^J 
+ 16 a*A*5 + 32 m^AA„(a)5 A^a^/3 + 24 rnX^& a 
-24m^A"*A/;„(a)a^ - 32a'* - 64m^(SA^a^ 



+ 16 A V 



96 rn?a^ 



-96 m^a^5 - 32 + 8 A* + 64 rnNm{a)\^a 
-32 m^Nm {a)\'afi + 16 a*A/m(a)5 A^/3 
+30AA™(q)/3 - 128mW,„(Q)a^/3 + 16AA™(a)<5Q^/3 
+32 m^AA™ (q)5 /3 + 32 a* Mm {a)S f3 
+64m^J\fmia)5a^P + 128 m^a^5 + 16 A/'„(a)5m^/3 
-112 - 30 5 /3 + 96 m*/3 + 96 Afm{a)m^P 
+ 112A/m(a)a^/3- 112 a^(5/3 + 128 m^a^/? 
-96 m''A/'„(Q)/3 + 2 AA™(a)(5 /3 + 2 A^/3 
+32m^5A^Q^/3 + 16 5A^Q^/3- 3A/'m(a)(5A^/3 
+32 m^A^Q^/3 + 32 A^/3 - 8 A/'m(a)A^m^/3 
-16A/'™(a)A'a^/3 - 32 a*AA™(a)A'/3 - 2AA™(a)A'/3 
+ 16 A^a^/3 + 32 a*A^/3 + 8 A^m^/3 + 14 5 + 32 a'*/? 
+32 a^S /3 - 32 q''aA™(q)/3 - SJ\fm{a)S a (3 A 
-64mW,„(a)5a^/3^ - 8 AA„(a)(5 A 
-l6Nra{a)5al3^ - 16AA„(a)5m^/3^ -32m^a^(5/3 
+32m2AA™(a)a^/32 - 32m''/3^ - 16AA„(a)/3 

+64a'^0^ + 16(5/3^ - 56 AA™(Q)m^/3^ 
-64A/'m(a)a^/3^ +64a^(5/3^ - 32 
+32m''AA„(a)/3^ - 2 A/'„(q)5/3^ + 40m^5 

-32 m*7V;„(a)5 - 32 Q'*AA™(a)5 + 56 

+40 + A^ - 30/3 + 48a*AA„(a)A^ 
+ 16/3^ - 8m''A''AA™(a) - 64m'' - A'' 
-Mm^S - 16a''A''AA™(a) + 9Qm^Nm{a)a 
+24 A*m^a^ + 96 m*(5 /3 + 8 m^<5 A^/3 - 96 rnS (3 
-32 771*5 56 m^(5 /3^ - 96 m^/3 + 32 a*A/'m(a) 
+64 77i'*AA„(q) + 25A^/3- 14AA„(a) . 



(D5) 



N 



|m 


=2, 8=0.2, a=1, ?i^=0.09 




p = 0.5589 







Figure El. The phase relationship between the perturbation of 
gas density fig and that of azimuthal magnetic field Z, with |77i| = 
2, (5 = 0.2, a = 1 and A'^ = 0.09. This branch goes across the 
horizontal j3 axis consistent with formula 11231 . 



parameters are exactly the same as those of Figure 1^ and 
only the y2 branch is shown here) to support our statement 
in the main text. In Figure El, it is clear that this y2 branch 
does indeed go across the horizontal /3 axis at about 0.5589 
as expected from formula p23^ . 



APPENDIX E: 

In section 4.3.2, we indicated that there exists a root /3cO for 
the 2/2 solution branch of the ^ig/Z, and in the case of Figure 
1231 the value of this /3cO is about 0.5589 by formula H123|l . 
However, one might worry that from Figure l^in which /3 
ranges from 1 to 5, the trend of variation of j/2 may not lead 
to such a /3co where y2 vanishes. Here, we provide the ^ig/Z 
versus (3 curve in which (3 ranges from 0.2 to 1.5 (all other 



